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1 . INTRODUCTION 


The primary purpose of this research effort was to develop 
personal computer based software for vibrational analysis. The 
software was developed to analytically determine the natural 
frequencies and mode shapes for the uncoupled lateral vibrations 
of the blade and counterweight assemblies used in a single bladed 
wind turbine. The uncoupled vibration analysis was performed in 
the flapwise and chordwise directions for static or non-rotating 
conditions of the counterweight assembly and in the flapwise 
direction for static conditions of the blade assembly. The 
effects of rotation on the uncoupled flapwise vibration of the 
blade and counterweight assemblies was evaluated for various rotor 
speeds up to 90 rpm. A secondary purpose of this research effort 
was to attempt to experimentally determine the natural frequencies 
for the blade and counterweight assemblies for static rotor 
conditions and to correlate these results with the analytical 
values . 


Section 2 of this report describes the theory used in 
calculating the natural frequencies and mode shapes. The theory 
is based upon a lumped mass formulation for the blade and 
counterweight assemblies. It is this theory which forms the basis 
for the development of the computer codes which ultimately 
generate the analytical results. Section 3 describes the computer 
codes developed and documents how the theory is implemented in the 
codes. The discussion of each code includes the input and output 
data structures used. The codes are designed to be as general as 
possible so that other designs can be readily analyzed. The input 
for the codes are generally interactive to facilitate usage. The 
actual code listings are included in the Appendix. Section 4 
describes the operation and use of the computer codes . Section 4 
also includes an example which demonstrates the application of the 
the computer software. Section 5 describes the actual results of 
the analysis of the MOD-O blade and counterweight assemblies at 
various operational speeds and includes a discussion of the 
analytical and experimental results. Section 6 contains general 
conclusions and recommendations. 



2. THEORY 


The following formulation concerns the uncoupled flapwise 
bending vibration of the blade and counterweight assemblies 
including the effects of rotation. Flapwise vibration, by 
definition, deals with bending in a plane perpendicular to the 
plane of rotation. Additional reading material on the 
determination of the vibratory characteristics of nonuniform beams 
and the effects of rotation on nonuniform rotating beams can be 
obtained in References (1) through (6). 


Figure 1 shows a discrete mass model of a nonuniform rotating 
cantilevered beam. The nonuniform beam has been modeled by a 
series of lumped masses connected by massless beam segments of 
various lengths. The flexural rigidity of each beam segment is 
allowed to vary along the length of the beam segment. Free body 
diagrams are illustrated in Figure 2 for mass n and the beam 
segments on either side of this mass. In Figure 2, 


l = the segment length 
m = the lumped mass 
M = the bending moment 
V = the shear force 

F = the centrifugal or inertia force 
due to rotation 


The deflection and slope of the beam at mass n can now be 
formulated in terms of the deflection and slope of the beam at the 
n-*l mass by letting 


d Sn = the deflection of a cantilevered beam at the 
free end due to a unit shear force applied 
at the free end 

d Mn = the deflection of a cantilevered beam at the 
>u free end due to a unit bending moment 
applied at the free end 




Figure 1 Lumped mass model of a nonuniform rotating beam. 






a Sn = t ^ le s ^°P e °f 3 cantilevered beam at the free 
1 end due to a unit shear force applied at the 
free end 

a Mn = t ^ le s *°P e a cantilevered beam at the free 
end due to a unit bending moment applied at 
the free end 


Expressions for these coefficients have been developed in detail 
in Reference (1). Assuming a linear variation in the flexural 
rigidity for a beam segment of length £ , these coefficients are 
given by 
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*n 

< EI >n " < EI) n-l 
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Where 


( EI ) n = the flexural rigidity of the beam 
at the location of mass n 

( EI ) n _l = the flexural rigidity of the beam 
at the location of the n-1 mass 


The centrifugal forces caused by rotation will remain parallel to 
the original undeformed axis of the beam as the beam vibrates 
perpendicular to the plane of rotation as illustrated in Figure 3. 
This effect will cause a stiffening of the beam, thus, increasing 
the natural frequencies of vibration. As seen in Figure 3, the 
centrifugal force on the right-hand end of the beam can be 
resolved into two components. Assuming small displacements and 
thus small slopes, one component is a tensile force acting on the 
beam of magnitude F E , and the other component is a downward 
shear force of magnitude F E 0 E . The deflection and slope of the 
beam segment at the location of mass n can then be determined in 
terms of the deflection and slope of the beam segment at the n-1 
mass location as illustrated in Figure 4. Thus, the deflection at 
the right-hand end of the beam segment is 


L R 

y = y . + 
-*n-l 


n 


0 R . + d M M L 
n-1 Mn n 


- d (vj; + fJ; ej*) 

Sn n n n 


( 2 ) 


and the corresponding slope is 


0 L = 0 R + a Mn mJ; - a (vj; + fJ; 0^) 
n n-1 Mn n Sn n n n 


(3) 


Taking moments about the left-hand end of the beam segment yields 


M L = M R + V L £ + F L 0 L £ 

n n-1 n n n n n 


Using the small slope assumption 


sin 0 L s tan 0^ = 0^ = 
n n n 


L R 
y n " y n-l 

*n 
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the bending moment in the beam segment at the location of mass n 
in the previous equation becomes 


M L = M R - + V L l + F L (y L - y R . ) 
n n-1 n n n ■'n J n-1 


(4) 


Since the beam segment shown in Figure 4 is massless, 
summation of the forces in the y-direction yields 


the 


V R - v L - F L e L + F L sin0 L = 0 

n-1 n n n n n 


Making use of the small deflection and 
assumption, this equation reduces to 


thus small slope 


V n =V nl 
n n-l 


(5) 


Using the free body diagram of mass n shown in Figure 2, Newton's 
law applied in the radial direction yields 


N 


F L = F R + m r ft 2 = S m. {r + x. 
n n n n / i o l 


)ft : 


( 6 ) 


i=n 


where, as shown in Figure 1, 


ft = the angular velocity of rotation 

r n = r G + x = the radial location of mass n 

r Q = the radius of the rigid hub 

N = the total number of lumped masses in 
the discrete model of the beam 


Applying Newton's second law in the radial direction 
body diagram of the massless beam segment of length 
Figure 2 yields 

N 

F R = F L = 3” m. (r + x . ) ft 2 
n-1 n / l o l 

i=n 


to the free 

& shown in 
n 


(7) 
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Using Equations 5 and 7, Equations 2, 3, and 4 become 


= 


■ n 


y R , + £ 0 R . + d M M L - d 
J n-1 n n-l Mn n Sn 




'+ F R e E ) 
n-l n 


( 8 ) 


0 L = 0 R i + a M M E - a q (V R + F R 0 B ) 
n n-l Mn n Sn n-l n-l n 


(9) 


m l = m r , 
n n-l 


+ *n V n-1 + F n-1 (F n ' **-!> 


(10) 


Following the method of formulation as presented in Reference 
(1), Equations 3 through 10 are solved for the deflection, slope 
and moment to the left of mass n in terms of the deflection, slope 
and moment to the right of mass n-*l . Writing the resulting 
equations in matrix form while utilizing equation 5, yields 
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[{ VW d Sn } *n (a Mn d Sn " a Sn d Mn } F n-1 ] /DEN 

R 


B 22 " [1 - < d Mn - V M n )F n-l ]/DEN 
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B 23 = a Mn /DEN 


24 


33 



a Mn 

■ a Sn> - (a Mn 

d Sn “ a Sn 

)F J/DEN 
Mn n-1 


‘32 = 

B 12 

F R , 
n-1 



(12 cont 

[1 + 

a Sn 

f b _ 1 j/den 




‘34 = 

B 12 





given by 





1 + 

(a Sn 

- d M )F R . + 
Mn n-1 

(a Mn d Sn 

a. d M ) (F R . ) 2 
Sn Mn n-1 

(13) 


To include the remaining dynamic effects associated with the 
vibration of the beam, one considers the free body diagram of the 
lumped mass illustrated in Figure 2. For continuity of the beam 
it follows that 



(14) 



(15) 


and 




(16) 


Summing the forces in the vertical direction yields 

V F = m 'y = V L - V R (17) 

/ y n^n n n 1 ' 

Assuming that each mass is subjected to harmonic motion in the 
vertical direction at the same frequency, 

y = A sin w t 
n 
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and 


y n = - o) 2 A n sin b) t = - w 2 y n (18) 


Substitution of y n from Equation 18 into Equation 17, the shear 
force acting to the right of mass n becomes 



Substituting for the deflection, slope, moment and shear to the 
left of mass n as given by Equation 11 into Equation 20 results in 
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Equation 21 can be written as 



( 22 ) 


where the coefficients of the A matrix are defined using Equations 
21, 12, 13 and 7. It should be noted that when the beam is not 
rotating. Equation 22 reduces to the equation for the vibration of 
a static or non- t rotating beam developed in Reference (1). Thus, 
Equation 22 can be used to investigate the chordwise and flapwise 
vibration of non-rotating beams by setting the angular velocity of 
the rotor to zero. If the beam is rotating, Equation 22 in its 
present form can only be used to investigate flapwise vibration. 


The determination of the natural frequencies follows the 
formulation presented in Reference (1). For the discrete model of 
the beam illustrated in Figure 5, Equation 22 can be applied to 
each of the lumped masses to yield, for n = 1 



(23) 
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station 0 station 1 station 2 station 3 


station N 



Figure 5 N mass model of a nonuniform rotating beam. 



and, for n = 2 



and, for n = 3 



and finally, for n = N 



Concatenating Equations 23 through 26 yields 


(24) 


(25) 


(26) 


15 





Letting the product of all the A matrices equal the matrix U f 
Equation 27 can be rewritten as 




(28) 


Since the beam shown in Figure 5 is considered to be fixed at 
station 0 and free at station N, the boundary conditions can be 
expressed as 
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Substituting the boundary conditions expressed by Equation 29 into 
Equation 28 and expanding the resulting matrix yields 


y N = u n (o) + u n<°) + u iV M * + u i4 v f 


12 


13 o 14 o 


(30) 


0 N u 21 (o) + u 22 (o) + U 23 M o + U 24 V o 


(31) 


0 


u 3 i ( °) + u 3 2 (o) + U 33 M o + U 34 V o 


(32) 


and 


u 4l(°) + u 42 (o) + u 43 M* + u 44 V* 


(33) 


For a non-trivial solution of Equations 32 and 33, that is, for 
the bending moment and shear force to the right of station 0 to be 
any value other than zero, the determinant of the coefficients of 
these terms must be equal to zero. Thus, 


or 


u 33 U 34 

U 43 U 44 


0 


u 


33 U 44 


u 


34 


u 


43 


(34) 


Equation 34 is the characteristic equation or frequency equation 
for the system. The roots of this equation, that is, the values 
of the frequencies that satisfy this equation, are the desired 
natural frequencies of vibration. 

Once the natural frequencies have been found using Equation 34 
, the corresponding mode shapes can be determined. The mode 
shapes be can obtained from the deflection of each of the lumped 
masses. In order to determine the deflection of the masses, 
Equation 33 can be solved for the bending moment to the right of 
station 0 yielding 
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( 35 ) 




u 


43 


Substituting the bending moment given by Equation 35 into Equation 
30, the deflection at the free end of the beam becomes 



u 


13 u 44 


u 


V R + 
o 


u 


14 


V* 


43 


(36) 


If the mode shapes are normalized with respect to the free end of 
the beam by letting the deflection at the free end be equal to 
one, the shear force at the fixed end of the beam obtained from 
Equation 36 becomes 



u 


43 


u 


14 U 43 


- u 


13 U 44 


(37) 


Substituting the shear force at the fixed end defined by Equation 
37 into Equation 35, the bending moment at the fixed end of the 
beam can be expressed as 



- u 


44 


U 14 U 43 “ U 13 U 44 


(38) 


Concatenating Equation 22 starting at mass 1 and going to mass n, 
the deflection, slope, bending moment and shear force to the right 
of mass n can be expressed in terms of the same quantities at the 
fixed end of the beam, station 0, as 



(39) 


Since the deflection and slope at the fixed end of the beam are 
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zero. Equation 39 can be rewritten as 




( 40 ) 


The normalized displacement for mass n can then be determined 
from Equation 40 by concatenating the A matrices as indicated and 
using the shear force and bending moment at station 0 as defined 
by Equations 37 and 38. The mode shape for a specific natural 
frequency is then obtained by determining the normalized 
deflection of each of the lumped masses. 
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3. DESCRIPTION OF COMPUTER CODES 


There are five computer codes which have been written for thv 
vibrational analysis described in Section 2.. Two codes, BLDFLAP 
and CWTFLAP , generate the discrete or lumped mass models used to 
analyze the flapwise vibration of the blade and counterweight 
assemblies respectively. An additional code, CWTCHRD , generates 
the lumped mass models used to analyze the chordwise vibration of 
the counterweight assembly. The only differences between the 
flapwise and chordwise codes are the changes in area moments of 
inertia due to the change in axes about which bending occurs. The 
fourth code, BEAM, is a general code for the determination of 
natural frequencies and mode shapes for cantilevered beams which 
have been modeled as lumped mass configurations. The fifth code, 
MODPLT, is a graphics code which displays the mode shapes on a 
Tektronix storage display terminal. The operation, data 
requirements, and output for each of the five codes and any 
associated subprograms are described in the following subsections. 


The analysis codes, written in FORTRAN 77, were originally 
developed on a Digital Equipment Corporation PDP 11/34 
minicomputer. The codes were then down loaded to an ITT XTRA 
personal computer. The ITT XTRA, a fully IBM compatible personal 
computer, utilized the Microsoft DOS 2.11 operating system and the 
Microsoft FORTRAN 3.20 compiler. After down loading, minor 
modifications were made to resolve FORTRAN incompatibilities and 
the codes were then compiled and linked. 


3.1 BLDFLAP 


The BLDFLAP computer code permits the user to input the number 
of lumped mannas in the various sections of the blade assembly. 
This assembly consists of the spool piece (CF 764254), shorty 40 
blade, transition section and 6 foot pitchahle tip. The program 
computes the distance between each of the lumped masses, the 
flexural rigidity at each lumped mass location and the weight of 
each lumped mass. Flexural rigidity is the product of the modulus 
of elasticity and area moment of inertia. The area moments of 
inertia are computed for bending in the flapwise direction. 
BLDFLAP then writes the title of the problem, the total number of 
lumped masses, the data for each lumped mass as described above, 
the flexural rigidity at the base of the spool piece arid the 
radius of the rigid hub to a data file whose name is " VIBIN .DAT" . 
The data which describes the blade geometry is defined in the code 
and was supplied by NASA (7), (8), (9). 
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3.2 CWTFLAP 


The CWTFLAP computer code permits the user to input the number 
of 1 urn pod inn.*: nos in thn various sections of the counterweight 
assembly. This assembly consists of the spool piece (CF 764254), 
the tapered stool, spar extension (CF 760549 ) and the counterweight 
(CF 764554). The computer code calculates the distance between 
each of the lumped mosses, the weight of each lumped mass and the 
flexural rigidity at each lumped mass location. The area moments 
of inertia used j.n the flexural rigidity calculations are computed 
for bending in the flapwise direction. CWTFLAP then writes the 
problem title, the total number of lumped mosses, the data for 
each lumped mass as previously described, the flexural rigidity at 
the base of the spool piece and the radius of the rigid hub to a 
data file whose name is " VIB1N .DAT" . The data which describes the 
counterweight geometry is defined in the code and was taken from 
the drawings as indicated. 


3.3 CWTCHRD 


The CWTCHRD computer code is identical to the CWTFLAP computer 
code described in Section 3.2 with the exception being that the 
area moments of inertia are taken assuming bending in the 
chordwise direction rather than in the flapwise direction. The 
only real difference is the moment of inertia of the mass at the 
free end of the counterweight assembly. The moment of inertia of 
this mass is not identical in the flapwise and chordwise 
direc tions . 


3.4 BEAM 


The BEAM code is the main vibration analysis code. It reads 
the lumped mass model data from the "VIBIN.DAT" data file and 
calculates the natural frequencies and the mode shapes. In 
addition to reading in the model data from the "VIBIN.DAT" data 
file, BEAM permits the user to enter the the range of values to be 
searched for natural frequencies. For flapwise vibration, the 
user may also specify the rotor speeds to be considered. The mode 
shape calculations are optional and thus must be requested if 
desired. The results of the frequency calculations are stored in 
a file called " VIBOUT .DAT"and the mode shape daLa are stored in a 
file entitled "MODKS.DAT". The BEAM main program is quite short 
with the bulk of the input, output and calculations being 
performed in eight subprograms. The eight subprograms are: DET , 
FTT.T.A , FORCE, TNPUT , MODES, MULT, NFREQ and OUTPUT. The BEAM main 
program and the eight subprograms are described in the following 
subsections below. 



3.4.1 BEAM 


The BEAM main program initially calls the INPUT subroutine 
which performs all the input functions. BEAM then enters a loop 
in which all natural frequencies and desired mode shapes are 
calculated and output for each of the requested rotor speeds. For 
chordwise vibration, only a static rotor condition is presently 
permitted. The natural frequencies are calculated by calling the 
NFREQ subroutine and are output to the "VIBOUT.DAT" data file by 
calling OUTPUT subroutine. The mode shapes are calculated and 
output to a file called "MODES. DAT” by calling the MODES 
subroutine . 


3.4.2 DET 


The DET function subprogram evaluates the determinant of the 
characteristic equation for the system. This determinant reflects 
the boundary conditions of the problem and is a function of the 
frequency. When the determinant is zero the frequency is a 

natural frequency. The actual determinant calculation is outlined 
in the theory section and is performed by first setting the U 
matrix to the identity matrix. Then the A matrix for each lumped 
mass is concatenated or multiplied with the existing U matrix 
This concatenation results in the U matrix for the entire beam. 
After the A matrices for all masses have been concatenated, the 
determinant, DET, is computed using the equation 

DET=U (3,3 )*U(4,4)-U(3,4)*U{4,3). 

The values of the frequency at which DET becomes zero are the 
natural frequencies of a cantilevered beam. 


3.4.3 FILLA 


The FILLA subroutine fills or defines the A matrix which 
relates the deflection, slope, moment and shear at mass n-1 to 
those quantities at mass n. The A matrices are functions of the 
frequency, centrifugal forces and geometric and material 
The exact form of the A matrix is given in the theory 
FILLA is called by DET prior to calculating the 


properties . 
section . 


determinant of the characteristic equation for the system. 
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3.4.4 FORCE 


The FORCE subroutine calculates the centrifugal force acting 
to the left of each mass for flapwise vibration of the beam. 
These values are stored in a force array and are then used later 
by F1LLA. This subroutine is called during the initial steps of 
the natural frequency calculation subroutine, NFREQ . If the non- 
rotating condition of the blade is being considered, the rotor 
speed which is passed to FORCE is 2 ero which results in all 
centrifugal forces being calculated as zero. 


3.4.5 INPUT 


The INPUT subroutine performs all input for the BEAM program. 
Initially INPUT reads the problem title, the total number of 
lumped masses, the lumped mass data, the flexural rigidity at the 
rigid hub and the radius of the rigid hub. All this data is read 
from the "V/1BIN.DAT" data file which was written by one of the 
data preparation programs such as CWTFLAP, BLDFLAP or CWTCHRD. 
The lumped mass data consists of the distance between each lumped 
mass, the flexural rigidity of the beam at each lumped mass 
location and the weight of each lumped mass. For flapwise 
vibration, INPUT then permits the user to interactively enter the 
number of rotor speeds to be considered and the values of these 
rotor speeds. The user then enters the lower and upper limits of 
the frequency range to be searched for natural frequencies and the 
increment to be used in the initial search. Finally the user 
enters the number of fractional decimal places of accuracy and 
indicates if mode shape calculations are desired. 


3.4.6 MODES 


The MODES subroutine calculates the mode shapes if such 
calculations were requested in the INPUT subroutine. Mode shapes 
are calculated for each of the natural frequencies at the various 
rotor speeds. The relative displacements at each lumped mass, or 
mode shapes, are output to a file call "MODES. DAT" for subsequent 
printing or plotting using the MODPLT program. The actual 
displacements for each natural frequency are calculated by again 
concatenating the A matrices for all lumped masses into a matrix U 
and then calculating the shear, VO, and moment, MO, at the fixed 
end assuming a unit displacement at the free end. The formulas 
for these calculations are 

MO = 1 ./ (U(l ,3 ) — U ( 4,3 ) *U ( 1 ,4 )/U( 4 ,4 ) ) 


and 


VO=-U (4,3) *MO/U ( 4,4). 
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Once the shear and moment are calculated at the fixed end, the 
displacements at each lumped mass, Y(N), can be calculated by 
simply concatenating A matrices up to the lumped mass whose 
displacement is desired and then apply the formula 

Y(N)=U(1 ,3 )*MO+U(l ,4)*VO. 


3.4.7 MULT 


The MULT subroutine is a simple matrix multiplication 
subroutine. It is used to concatenate the A matrices in the DET 
and MODES subprograms. 


3.4.8 NFHEQ 


The NFREQ subroutine computes the natural frequencies in the 
range of frequencies as specified in the INPUT subroutine. Before 
the subroutine can actually calculate the natural frequencies it 
first computes the centrifugal force at each lumped mass by 
calling the FORCE subroutine. Once this is done, the NFREQ 
subroutine then begins to search for natural frequencies using the 
method of interval halving. The natural frequencies are those 
frequencies at which the characteristic determinant of the system 
is zero. This interval halving method calls on the DET function 
to evaluate the determinant as described earlier. When a natural 
frequency is found to the precision requested in the INPUT 
subroutine, it is stored in an array and the search continues for 
the next possible frequency in the range of the search. The 
natural frequencies are ultimately written by the OUTPUT 
subroutine into the file " V1B0UT .DAT" . The NFREQ program, 
however, also prints the natural frequencies at the terminal so 
the user can see the progress the program is making. 


3.4.9 OUTPUT 


The OUTPUT subroutine outputs the input data and nalural 
frequencies to the file "VIBOUT.DAT" for subsequent printing. The 
output is formatted and labelled to facilitate comprehension. The 
input data is echoed on the printed output so that the listing of 
"VIBOUT.DAT" provides one complete summary of the analysis 
including results and the original discrete model. 
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3.5 MODPLT 


MODPLT is a stand alone program which plots the mode shapes on 
a Tektronix 4010 or 4014 display terminal. The program reads the 
"MODES.DAT" data file and permits the user to generate a quick 
sketch, to scale the results in order to create a pleasing plot, 
or to generate a finished labelled plot. The program was written 
with all the graphics subroutine calls in one subroutine called 
DRAW. This will hopefully minimize the difficulty in moving this 
program to another machine. The only other subroutine used in 
MODPLT is a subroutine call MINPlfT which actually performs the 
input from "MODES.DAT". 
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4. OPERATION AND USE OF COMPUTER CODE 


This section illustrates the use of the programs discussed in 
Section 3 to generate a discrete model, to perform a vibration 
analysis and to plot the corresponding mode shapes. Section 4.1 
1 1 1 us h rates the use of the BLDFJ.AP program to generate a 107 
lumped mass model of the blade assembly for flapwise vibration. 
Included in this section is a hard-copy of the interaction between 
the computer and the user taken directly from the computer screen. 
As in all sections to follow where computer- user interaction is 
shown, the specific user response has been underlined for the 
reader's convenience . Section 4.2 contains a listing of the VIBIN 
data file generated by the BLDFLAP program. Section 4.3 
illustrates the use of the BEAM program to perform a flapwise 
vibration analysis of the blade assembly for rotor speeds of 0, 20 
and 40 revolutions per minute. Included in this section is a 
hard-copy of the computer- user interaction. The results of this 
analysis are the first three natural frequencies of the blade 
assembly and the corresponding mode shapes based on the discrete 
model developed in Section 4.1. Section 4.4 contains a listing of 
the V1BOUT data file generated by the BEAM program . This data 
file includes a complete listing of all input data and the values 
of the natural frequencies as found by the BEAM program. Section 
4.5 contains a listing of the MODES data file generated by the 
BEAM program during the vibration analysis. This data file 
includes all pertinent information relative to the corresponding 
mode shapes for the various natural frequencies of vibration. 
Section 4.6 illustrates the use of the MODPLT program to plot the 
mode shapes on a Tektronix 4010 or 4014 display terminal . Included 
in this section is a hard-copy of the computer-user interaction 
and copies of the three mode shapes for each of the three rotor 
speeds used in the analysis. 



4.1 EXAMPLE OF MODEL CREATION 


The following example illustrates the use of the BLDFLAP 
program which generates the discrete model of the blade assembly. 
This model can be used for a flapwise vibration analysis performed 
by the BEAM program. The program is executed by entering BLDFLAP. 
The program then queries the user for the radius of the rigid hub, 
the number of divisions for each section of the blade and the 
title to be assigned to the problem. The values shown in the 
example are for illustrative purposes only and may be altered to 
suit the user's requirements. The total weight and the location 
of the center of gravity of the discrete model are displayed for 
the user 's convenience as well as the individual weight and 
location of the center gravity for each section of assembly. 

This program is typical of the operation of the other model 
generation programs , CWl'FLAP and CVTCHRD . Thus, illustrations for 
these two programs was deemed unnecessary. The input for all 
three model generation programs is totally interactive and the 
output is the V1BJN data file. This output file is denoted as 
VIBIN because it is used as input to the analysis code, BEAM. The 
interactive input follows and the output data file is shown in 
Section 4.2. 
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B> BLDFLAP 


THIS PROGRAM GENERATES THE DATA FOR THE FLAPWISE 
VIBRATION OF THE BLADE ASSEMBLY 

ENTER THE RADIUS OF THE RIGID HUB (INCHES) : ££. O 

SECTION 1-SPOOL PIECE FLANGE L=£. £5” W=108.74 LB 
SECTION £— SPOOL PIECE BODY L=13. 75" W=£07. 96 LB 
SECTION 3-SPOOL PIECE FLANGE L=1.75" W=S£. 78 LB 
SECTION 4-BLADE SECTION L=441.£5" W=££50. 00 LB 
SECTION 5-TRANSITION & TIP SECTION L=115. 58" W=960. 00 LB 

00 

THE MAXIMUM TOTAL NUMBER OF DIVISIONS YOU 
ARE ALLOWED IN ALL 5 SECTIONS IS 699. 


ENTER THE NUMBER OF DIVISIONS FOR SECTION 1 : £ 


ENTER THE NUMBER OF DIVISIONS FOR SECTION £ : 3 


ENTER THE NUMBER OF DIVISIONS FOR SECTION 3 : £ 


ENTER THE NUMBER OF DIVISIONS FOR SECTION 4 : 80 


ENTER THE NUMBER OF DIVISIONS FOR SECTION 5 : £0 



FLAPWISE VIBRATION OF THE BLADE 


ENTER A TITLE STATEMENT (60 CHAR. MAX. ) : 


TOTAL WEIGHT 
(LBS) 
3569. 472 


SPOOL WEIGHT 
(LBS) 

399. 472 


BLADE WEIGHT 
(LBS) 
2250. 000 

N> 

<o 

TRANS. & TIP 
(LBS) 

920. 000 


CG FROM HUB END 
(INCH) 

266. 502 


CG FROM HUB END 
(INCH) 

8.553 


CG FROM SPOOL END 
(INCH) 

200. 300 


WEIGHT CG FROM BLADE END 

(INCH) 

38. 004 


TOTAL MOMENT ABOUT CENTERLINE OF HUB 
(LBS-INCH) 

1029800. 5 

Stop - Program terminated. 



4.2 EXAMPLE OF VIBIN DATA FILE 


The following is a listing of the VIBTN data file generated by 
the BT.DFLAP program for the example problem illustrated in Section 
4.1. This data file contains a complete description of the 
problem including: the problem title, the number of lumped masses, 
the section lengths, the section rigidities and the lumped 
weights. The VTBTN data file is an intermediate file and all the 
information in this file is contained in the analysis code output 
data file which is formatted with titles suitable for printing. 
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cu <r o- oj 


FLAPWISE VIBRATION OF THE BLADE 


107 

. 562500 
1. 125000 
. B54167 
. 583333 

• 5Bi^^fbi3 

. 729167 
. 875000 
3. 195313 
5. 527500 
5. 503750 
5. 515625 
5.515625 
5. 515625 
5. 006254 
5. 944996 
5. 515625 
5. 515625 
5. 515059 
5. 516191 
5. 515625 
5. 515625 
5. 515625 
5. 348701 
5. 602549 
5. 515625 
5. 515625 
5. 515347 
5. 515903 
5. 515625 
5. 515625 
5. 515625 
5. 478006 
5. 553164 
5. 515625 
5. 515625 
5. 515625 


229986015274. 6 
229986015274. 6 
63188041992. 2 
63188041992.2 

63188041992. 2 
226554261922. 3 
226554261922. 3 

28147549273. 9 
27441130208. 7 

26737746369. 3 
26032844917. O 
25327943464. 7 
24623042012. 4 

23973014492. 5 

232 1 3239 1 07 . 9 

22508337655. 6 

21803436203. 3 

21098607134. 7 

20504289062. 5 

19991335937. 5 

19478382812. 5 

18965429687. 5 

18468000460. 9 

17939523437. 5 

17426570312. 5 

16913617187. 5 

16400689925. 2 

15887710937. 5 

15374757812. 5 

14861804687.5 

14348851562.5 

13839389563. 3 

13322945312. 5 

12809992187. 5 

12297039062. 5 

11784085937. 5 


54. 

54. 

69. 

69. 

69. 

41. 

41. 

61. 

62. 

64. 

64. 

64. 

64. 

45. 

34. 

34. 

34. 

34. 

34. 

34. 

34 

34 

30 

26. 

26 

26 

26 

26 

26 

26 

26 

25 

25 

25 

25 

25 


367651 

367651 

319249 

319249 

319249 

389338 

389338 

952576 

264924 

157750 

157750 

157750 

157750 

884094 

154956 

154956 

154956 

147363 

055675 

055675 

055675 

055675 

065675 

719894 

719894 

719894 

717169 

479413 

479413 

479413 

479413 

805525 

071825 

071825 

071825 

071825 



u> 

N> 


5. 504605 

11377504100. 4 

23 . 339200 

5. 526645 

11101171875.0 

23 . 240638 

5. 515625 

10825390625. 0 

23 . 240638 

5. 515625 

10549609375. 0 

23. 240638 

5. 427723 

10278223203. 5 

22. 042525 

5. 603527 

9998046875. 0 

20. 343831 

5. 515625 

9722265625. 0 

20. 343831 

5. 515625 

9446484375. 0 

20. 343831 

5. 515625 

9170703125. 0 

20. 343831 

5. 482032 

8896601504. 8 

18. 454506 

5. 549218 

8619140625. 0 

18. 199356 

5. 515625 

8343359375. 0 

18. 199356 

5. 515625 

8067578125. 0 

18. 199356 

5. 416897 

7796733265. 7 

17. 251488 

5. 614353 

7516015625. 0 

15. 479050 

5. 515625 

7240234375. 0 

15. 479050 

5. 515625 

6964453125. 0 

15. 479050 

5. 515625 

6688671875. 0 

15. 479050 

5. 545001 

6475198116. 6 

16. 429450 

5. 486249 

6398390625. 0 

16. 643950 

5. 515625 

6321171875. 0 

16. 643950 

5. 515625 

6243953125. 0 

16. 643950 

5. 578573 

6165853102. 5 

17. 191450 

5. 452677 

6089515625. 0 

18. 576625 

5. 515625 

6012296875. 0 

18. 576625 

5. 515625 

5935078125. 0 

18. 576625 

5. 515625 

5857859375. 0 

18. 576625 

5. 412902 

5782078741. 6 

16. 161813 

5. 618348 

5703421875. 0 

15. 359913 

5. 515625 

5626203125. 0 

15. 359913 

5. 515625 

5548984375. 0 

15. 359913 

5. 621312 

5470286008. 8 

16. 151569 

5. 409938 

5394546875. 0 

18. 989194 

5. 515625 

5317328125. 0 

18. 989194 

5. 515625 

5240109375.0 

18. 989194 

5. 515625 

5162890625. 0 

18. 989194 

5. 597277 

5078144094. 5 

20. 963006 

5. 433973 

4970673671. 2 

21. 868350 



5. 515625 

4861588380 . 7 

21. 868350 

5. 515625 

4752503090 . 2 

21. 868350 

5. 582987 

4642085544 . 3 

22. 517700 

5. 448263 

4534332509 . 3 

26. 113175 

5. 515625 

4425247218 . 8 

26. 113175 

5. 515625 

4316161928 . 3 

26. 113175 

5. 515625 

4207076637 . 8 

26. 113175 

5. 809453 

4092180165 . 8 

36. 302150 

5. 221797 

3988906056 . 9 

42. 536500 

5. 515625 

3879820766. 4 

42. 536500 

5. 515625 

3770735475. 9 

42. 536500 

5. 765264 

3656712941. 5 

47. 222667 

5. 265986 

3552564894. 9 

95. 898333 

5. 646812 

9133101005. 6 

99. 138314 

5. 778000 

8607303016. 7 

99. 138314 

5. 778000 

8384209507. 0 

99. 138314 

5. 596266 

8253244818. 2 

81. 383133 

5. 959734 

7710909050 . 0 

74. 189518 

5. 778000 

7185111061 . 1 

74. 189518 

5. 778000 

6659313072 . 2 

74. 189518 

5. 129401 

1 952000000 . 0 

54. 351447 

6. 426598 

1511815488 . 7 

22. 649760 

5. 778000 

1471943302 . 3 

22. 649760 

5. 778000 

1432522408 . 7 

22. 649760 

5. 766651 

1393284319 . 7 

22. 555680 

5. 789349 

1353337810 . 5 

21. 031920 

5. 778000 

1313651211. 2 

21.031920 

5. 778000 

1274360812 . 0 

21. 031920 

5. 778000 

1234733213 . 1 

21. 031920 

5. 815978 

1194602963 . 2 

22. 114214 

5. 740022 

1155208014 . 6 

22 . 51 1688 

5. 778000 

707573008 . 8 

22 . 51 1688 

5. 778000 

235857802 . 5 

22 . 51 1688 


229986015274 . 6 
22. 000 



4.3 EXAMPLE OF VIBRATION ANALYSIS 


The following illustrates the use of the BEAM program to 
perform a vibration analysis of the discrete model creaked by the 
BLDFLAP program illustrated in Section 4.1. The program is 
executed by entering BEAM. The program then queries the user for 
the type of analysis to be performed (chordwise or flapwise 
vibration), the rotor speeds to be considered {provided a flapwise 
analysis is requested), the frequency range and search interval, 
the precision and whether or not the mode shapes are to be 
determined. For the user's convenience, the program displays the 
natural frequencies at the various rotor speeds ns they are 
calculated . 



B> BEAM 

TYPE 1 IF THE ANALYSIS IS FOR CHORDWISE VIBRATION OR 
TYPE £ IF THE ANALYSIS IS FOR FOR FLAPWISE VIBRATION : £ 

TYPE IN THE TOTAL NUMBER OF ROTATIONAL 


SPEEDS 

TO BE l 

CONSIDERED (MAX. 

10) 

: _3 


TYPE 

IN 

ROTOR 

SPEED 

( REV/MIN) 

NO. 

1 : 

0. 0 

TYPE 

IN 

ROTOR 

SPEED 

( REV/MIN) 

NO. 

£ : 

O 

■ 

O 

OJ 

TYPE 

IN 

ROTOR 

SPEED 

( REV/MIN) 

NO. 

3 : 

o 

• 

o 


TYPE IN THE LOWER LIMIT OF THE FREQUENCY RANGE 

YOU WISH SEARCHED FOR NATURAL FREQUENCIES <RAD/SEC> : 0. O 

TYPE IN THE UPPER LIMIT OF THE FREQUENCY RANGE 

YOU WISH SEARCHED FOR NATURAL FREQUENCIES (RAD/SEC) : ISO. 

TYPE IN THE FREQUENCY INCREMENT TO BE USED IN 

THE INITIAL SEARCH FOR NATURAL FREQUENCIES <RAD/SEC> : 10. 

TYPE IN THE NUMBER OF DIGITS OF ACCURACY TO THE 
RIGHT OF THE DECIMAL POINT : 3 


DO YOU WISH TO CALCULATE MODE SHAPES (Y/N) ? Y 



** THE PROGRAM IS NOW DETERMINING THE NATURAL 
FREQUENCIES OF VIBRATION ** 

NATURAL FREQUENCIES AT A ROTOR SPEED OF .0 RPM 

************************************************** 

W = 10.341 RAD/SEC 

67. 686 RAD/SEC 

158. 607 RAD/SEC 

** THE PROGRAM IS CALCULATING THE MODE SHAPES ** 




** THE PROGRAM IS NOW DETERMINING THE NATURAL 
FREQUENCIES OF VIBRATION ** 


NATURAL FREQUENCIES AT A ROTOR SPEED OF £0.0 RPM 
************************************************** 



10. 650 
68. 066 
158. 909 


RAD/SEC 

RAD/SEC 

RAD/SEC 


** THE PROGRAM IS CALCULATING THE MODE SHAPES ** 



** THE PROGRAM IS NOW DETERMINING THE NATURAL 
FREQUENCIES OF VIBRATION ** 

NATURAL FREQUENCIES AT A ROTOR SPEED OF 40.0 RPM 
•************************************************** 

W = 11.5SO RAD/SEC 

W = 69. 194 RAD/SEC 

W = 159. BOS RAD/SEC 

** THE PROGRAM IS CALCULATING THE MODE SHAPES ** 
Stop - Program terminated. 



4.4 EXAMPLE OP VIBOUT DATA PILE 


The following is a listing of the VIBOUT data file generated 
by the BEAM program for the example problem illustrated in the 
previous sections. This data file contains a complete description 
of the problem including: the problem title, the number of lumped 
masses, the section lengths, the section rigidities, the lumped 
weights, the section rigidity at the base of the model, the radius 
of the rigid hub, the frequency range and search interval, and the 
natural frequencies at the various rotor speeds considered. The 
VIBOUT data file can be spooled to a printer in order to provide 
the user with a hard-copy of all input and output relative to the 
vibration analysis performed. 



FLAPWISE VIBRATION OF THE BLADE 


THE NUMBER OF LUMPED MASSES » 107 


SECTION LENGTH 
(INCH) 

. 562500 
1. 125000 

© 2.854167 

4. 583333 
4. 583333 
2.729167 
. 875000 
3. 195313 
5. 527500 
5. 503750 


SECTION RIGIDITY 

LUMPED WEIGHT 

(LBS-IN*IN> 

(LBS) 

229986015274. 6 

54. 367651 

229986015274. 6 

54. 367651 

63188041992.2 

69. 319249 

63188041992. 2 

69. 319249 

63188041992. 2 

69. 319249 

226554261922. 3 

41. 389338 

226554261922. 3 

41. 389338 

28147549273.9 

61. 952576 

27441130208.7 

62. 264924 

26737746369. 3 

64. 157750 



5. 515625 


5. 515625 
5. 515625 
5. 086254 
5. 944996 
5. 515625 
5. 515625 
5. 515059 
5. 516191 
5. 515625 
5. 515625 
5. 515625 
5. 348701 
5. 682549 
5. 515625 
5. 515625 


5. 515347 


26032844917. 0 

64. 157750 

25327943464. 7 

64. 157750 

24623042012 . 4 

64. 157750 

23973014492.5 

45. 884094 

23213239107 . 9 

34. 154956 

22508337655. 6 

34. 154956 

21803436203 . 3 

34. 154956 

21098607134. 7 

34. 147363 

20504289062. 5 

34. 055675 

19991335937. 5 

34. 055675 

19478382812. 5 

34. 055675 

18965429687. 5 

34. 055675 

18468000460. 9 

30. 065675 


17939523437.5 26.719894 

17426570312.5 26.719894 

16913617187.5 26.719894 

16400689925.2 26.717169 



5. 315903 

15BB77 10937. 5 

26. 479413 

5. 515625 

15374757812. 5 

26. 479413 

5. 515625 

14861804687. 5 

26. 479413 

5. 515625 

14348851562. 5 

26. 479413 

5. 47B086 

13839389563. 3 

25. 805525 

5. 5531 64 

13322945312.5 

25. 071825 

5. 515625 

12809992187. 5 

25. 071825 

5. 515625 

12297039062.5 

25. 071825 

5.515625 

11784085937.5 

25. 071825 

5. 504605 

11377504100. 4 

23. 339200 

5. 526645 

11101171875. 0 

23. 240638 

5. 515625 

10825390625. 0 

23. 240638 

5. 515625 

10549609375.0 

23. 240638 

5. 427723 

10278223203. 5 

22. 042525 

5. 603527 

9998046875. 0 

20. 343831 

5. 515625 

9722265625. 0 

20. 343831 

5.515625 

9446484375. 0 

20. 343831 



5. 515625 
5. 482032 
5. 549218 

5. 515625 
5. 515625 
5. 416897 
5. 614353 
5. 515625 

5.515625 

u> 

5. 515625 
5. 545001 
5. 486249 
5. 515625 

5. 515625 
5. 578573 
5. 452677 
5. 515625 

5. 515625 


9170703125. 0 

20. 343831 

8896601504. 8 

18. 454506 

8619140625. 0 

18. 199356 

8343359375. 0 

18. 199356 

8067578125. 0 

18. 199356 

7796733265. 7 

17. 251488 

7516015625. 0 

15. 479050 

7240234375. 0 

15. 479050 

6964453125. O 

15. 479050 

6688671875. 0 

15. 479050 

6475198116. 6 

16. 429450 

6398390625. 0 

16. 643950 

6321171875. 0 

16. 643950 

6243953125. 0 

16. 643950 

6165853102. 5 

17. 191450 

6089515625.0 

18. 576625 

6012296875.0 

18. 576625 

5935078125. 0 

18. 576625 



*- 

■> 


5. 515625 
5. 412902 
5. 618348 

5.515625 

5. 515625 
5. 621312 
5. 409938 

5.515625 

5. 515625 

5. 515625 
5. 597277 
5. 433973 

5. 515625 

5.515625 
5. 582987 
5. 448263 


5857859375. 0 

18. 

5782078741. 6 

16. 

5703421875.0 

in 

5626203125.0 

15. 

5548984375. 0 

15. 

5470286008. 8 

16. 

5394546B75. 0 

18. 

5317328125.0 

18. 

5240109375.0 

18. 

5162890625.0 

18. 

5078144094.5 

20. 

4970673671.2 

21. 

4861588380.7 

21. 

4752503090. 2 

21. 

4642085544. 3 

22. 

4534332509. 3 

26. 

4425247218. 8 

26. 


576625 

161813 

359913 

359913 

359913 

151569 

989194 

989194 

989194 

989194 

963006 

868350 

868350 

868350 

517700 

113175 

113175 


5. 515625 



5. 515625 


5. 515625 
5. 809453 
5. 221797 

5. 515625 

5. 515625 
5. 765264 
5. 265986 
5. 646812 

5. 778000 

5. 778000 
5. 596266 
5. 959734 

5. 778000 

5. 778000 
5. 129401 


6. 426598 


4316161928. 3 

26. 113175 

4207076637. 8 

£ 6. 113175 

4092180165.8 

36. 302150 

3988906056. 9 

42. 536500 

3879820766. 4 

42. 536500 

3770735475. 9 

42. 536500 

3656712941. 5 

47 . 222667 

3552564894 . 9 

95. 898333 

9133101005 . 6 

99. 138314 

8607303016 . 7 

99. 138314 

8384209507. 0 

99. 138314 

8253244818 . £ 

81. 383 1 33 

7710909050. 0 

74. 189518 

7185111061 . 1 

74. 189518 

6659313072 . £ 

74. 189518 

1952000000. 0 

54. 351447 


1511815488 . 7 


22. 649760 



5. 778000 

1471943302. 3 

22. 649760 

5. 778000 

1432522408.7 

22. 649760 

5. 766651 

1393284319. 7 

22. 555680 

5. 789349 

1353337810. 5 

21. 031920 

5. 778000 

131365121 1. 2 

21. 031920 

5. 778000 

1274360812. O 

21. 031920 

5. 778000 

1234733213. 1 

21. 031920 

5. 815978 

1194602963. 2 

22. 114214 

5. 740022 

1155208014. 6 

22. 511688 

5. 778000 

707573008. 8 

22.511688 

5. 778000 

235857802. 5 

22.511688 


THE SECTION RISIDITY (AT X=0. O) = 229986015274. 6 (LBS-IN*IN> 

THE RADIUS OF THE RIGID HUB = 22.000 (IN) 

THE FREQUENCY RANGE SEARCHED FOR NATURAL FREQUENCIES STARTING AT: 

00 AND ENDING AT: 180.00000 (RAD/SEC) 

THE INITIAL FREQUENCY INCREMENT USED IN THE SEARCH = 10.00000 


. 000 

(RAD/SEC) 



•t* 


NATURAL FREQUENCIES AT A ROTOR SPEED OF .0 RPM 

********************************-****************** 

(RAD/SEC) <HZ) 


1 10.341 1.646 

£ 67.686 10.773 

3 158.607 £5. £43 


NATURAL FREQUENCIES AT A ROTOR SPEED OF £0.0 RPM 
***********************-*************************** 

(RAD/SEC) <HZ) 


1 10.650 1.695 

£ 68.066 10.833 


3 


158. 909 


£5. £91 



J* 

00 


NATURAL FREQUENCIES AT A ROTOR SPEED OF 40.0 RPM 
************************************************** 

(RAD/SEC) (HZ) 

1 11.520 

£ 69. 194 

3 159.808 


1. 834 
11. 013 
£5. 434 



4.5 EXAMPLE OF MODES DATA FILE 


The following is a listing of the MODES data file generated by 
the BEAM program for the example illustrated in Sections 4.1 and 
4.3. This data file contains a complete description of the 
problem including: the problem title, the total number of lumped 
masses in the beam model, the axial location of each mass, the 
number of natural frequencies determined for each of the rotor 
speeds, the natural frequencies of vibration and the normalized 
displacements for each of the lumped masses at the natural 
frequencies. The MODES data file is used by the MODPLT program in 
order to obtain plots of the mode shapes. 
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THE PROBLEM TITLE 
***************** 


FLAPWISE VIBRATION OF THE BLADE 


THE TOTAL NUMBER OF MASSES 
************************** 

107 


ui 

o 


THE AXIAL LOCATION OF MASS N (INCHES) 
************************************* 


. 56250 
£ 6. 03531 
70. 14844 
114. £7344 
158. 39844 
£ 0£. 43554 
£ 46. 64844 
£ 90. 77344 
334. 89844 
379. 10509 
4 £ 3. 14844 
467. 66700 
513. 89100 
560. 11500 


1. 68750 
31. 53906 
75. 66350 
119. 78906 
163. 91406 
£ 08. 03906 
£ 5£. 06534 
£ 96. £8906 
340. 41406 
384. 53906 
4 £ 8. 95789 
473. 44500 
519. 66900 
565. 89300 


4. 54167 
37. 05469 
81. 17969 
1 £ 5. 30441 
169. 42969 
£ 13. 55469 
£ 57. 67969 
301. 86764 
345. 92969 
390. 05469 
434. 17969 
479. 04127 
5 £ 5. 43565 
571. 67100 


9 . 1£500 
4£. 57031 
86. 6953 1 
130. 82031 
174. 94531 
£ 19. 07031 
£ 63. 19531 
307. 32031 
351. 55100 
395. 57031 
439. 69531 
485. 00100 
53 1 . ££500 


13. 70833 
48. 08594 
9£. £1094 
136. 33594 
180. 44992 
£ 24. 58594 
£ 68. 71094 
31£. 83594 
356. 96094 
401. 15330 
445. £1094 
490. 77900 
537. 00300 


16. 43750 
53. 17219 
97. 72656 
141. 85156 
185. 97656 
£ 30. 06797 
£ 74. ££656 
318. 35156 
362. 47656 
406. 60156 
450. 97620 
496. 55700 
542. 78100 


17. 31250 
59. 11719 
103. 07526 
147. 36719 
191. 49219 
£ 35. 61719 
279. 77156 
323. 86719 
367. 99219 
41£. 1 1719 
456. £4219 
501. 68640 
548. 55900 


£ 0. 50781 
64. 63281 
108. 75781 
152. 84527 
197. 00781 
£ 41. 13281 
£ 85. £5781 
329. £8009 
373. 50781 
417. 63281 
461. 88900 
508. 11300 
554. 37498 



THE NUMBER OF NATURAL FREQUENCIES AT THE LISTED ROTOR SPEED (RPM) 
******************************************************************** 

3 .0 


THE NATURAL FREQUENCY OF VIBRATION <RAD/SEC> 
******************************************** 

10.34122 


THE NORMALIZED DISPLACEMENTS AT THE LUMPED MASSES 
**********■*****■**•********#**■*■**«•■**■*■*■*•*#••*•*•#*•***■#** 


00000 

. 00000 

.00001 

. 00005 

. 00014 

. 00022 

. 00024 

. 00034 

00061 

. 00102 

.00160 

. 00233 

. 00323 

. 00420 

.00551 

. 00690 

00847 

. 01021 

. 01213 

. 01423 

. 01651 

. 01897 

. 02154 

. 02446 

02749 

. 03072 

. 03414 

. 03776 

. 04158 

. 04560 

. 04984 

. 05425 

05894 

. 06383 

. 06893 

. 07427 

. 07983 

. 08565 

. 09170 

. 09799 

10442 

. 11130 

. 11832 

. 12558 

. 13310 

. 14081 

. 14888 

. 15714 

16567 

. 17429 

. 18349 

. 19279 

. 20235 

. 21219 

. 22234 

. 23266 

24331 

. 25421 

. 26550 

. 27679 

. 28845 

. 30035 

. 31249 

. 32462 

33745 

. 35026 

. 36329 

. 37678 

. 38996 

. 40360 

.41742 

. 43142 

44582 

. 45996 

. 47447 

.48915 

. 50416 

. 51895 

. 53407 

. 54932 

56469 

. 58101 

. 59579 

. 61151 

. 62732 

. 64394 

.65921 

. 67564 

69248 

. 70936 

. 72572 

. 74317 

. 76010 

. 77705 

.79211 

. 81 103 

82808 

. 84518 

. 86227 

. 87947 

. 89665 

. 91385 

.93107 

. 94841 

96553 

. 98276 

1. 00000 








THE NATURAL FREQUENCY OF VIBRATION (RAD/SEC) 
****************************-* -a********-# ■*•**•** 

67. 68574 


THE NORMALIZED DISPLACEMENTS AT THE LUMPED MASSES 
************************************************* 


00000 

00001 

00005 

00032 

-. 00090 

-. 00135 

-. 00151 

-. 00211 

00370 

00617 

00952 

-.01374 

-. 01883 

-. 02426 

-. 03152 

-. 03911 

04751 

05670 

06668 

07742 

08890 

-. 10109 

-. 11358 

-. 12752 

14170 

15649 

17185 

-. 18776 

-. 20419 

-. 22109 

-.23843 

-. 25605 

27428 

29271 

31142 

— . 33036 

-. 34944 

-. 36873 

-. 38805 

-. 40738 

42636 

44584 

46485 

-. 48362 

-.50210 

-. 52009 

-. 53786 

-. 55502 

57159 

58722 

60266 

61700 

-. 63043 

-. 64286 

-. 65424 

-. 66433 

67323 

68081 

68708 

69180 

-. 69508 

-. 69683 

-. 69698 

-. 69553 

69232 

68742 

68076 

— . 67212 

— . 66200 

-. 64984 

-. 63579 

-. 61983 

60166 

58210 

56029 

53649 

-. 51036 

-. 48288 

-. 45305 

-. 42120 

38735 

34953 

31365 

— . 27386 

-. 23214 

-. 18651 

-. 14308 

-. 09500 

04498 

. 00572 

. 05545 

. 10903 

. 16156 

. 21463 

. 26224 

. 32339 

38014 

. 43842 

. 49789 

. 55872 

. 62035 

. 68271 

. 74563 

. 80936 

87251 

. 93621 

1 . 00000 








THE NATURAL FREQUENCY OF VIBRATION (RAD/SEC) 
******************************************** 

15B. 60729 


THE NORMALIZED DISPLACEMENTS AT THE LUMPED MASSES 

♦t********###*#*#*###*******#******#***#*#*#**#*# 


00000 

. 00001 

. 00010 

. 00069 

. 00189 

. 00285 

. 00318 

. 00443 

00770 

. 01269 

. 01937 

. 02767 

. 03750 

. 04788 

. 06152 

. 07554 

09080 

. 10721 

. 12468 

. 14310 

. 16237 

. 18238 

. 20239 

. 22417 

24571 

. 26752 

. 28947 

. 31144 

. 33329 

. 35490 

. 37612 

. 39668 

41684 

. 43606 

. 45431 

. 47146 

. 48732 

. 50183 

. 51480 

. 52612 

53555 

. 54340 

. 54917 

. 55289 

. 55449 

. 55390 

. 55103 

. 54586 

53833 

. 52860 

. 51605 

. 50127 

. 48405 

. 46442 

. 44226 

.41802 

39145 

. 36280 

. 33187 

. 29992 

. 26605 

. 23082 

. 19444 

. 15784 

11916 

. 08075 

.04218 

. 00298 

03437 

07177 

10820 

-. 14336 

17742 

20861 

23807 

26498 

28927 

30980 

32705 

-. 34040 

34954 

35425 

35388 

34847 

33763 

32025 

29875 

-. 27055 

23877 

20438 

16861 

12795 

08603 

04167 

. 00002 

. 05928 

12099 

80596 

- 19013 
. 90276 

. 26563 
1 . OOOOO 

. 34703 

. 43296 

. 52265 

. 61525 

. 71056 



THE NUMBER OF NATURAL FREQUENCIES AT THE LISTED ROTOR SPEED (RPM) 
**************************************#*************«-************ 

3 20. O 


THE NATURAL FREQUENCY OF VIBRATION (RAD/SEC) 
******************************************** 

10. 649SB 


THE NORMALIZED DISPLACEMENTS AT THE LUMPED MASSES 
***#*************«■■*•****■***■#****■**■**■•#•**■#•*•**•*■*■***•** 


00000 

. 00000 

. 00001 

. 00005 

.00015 

. 00022 

. 00025 

. 00035 

00062 

.00104 

. 00163 

. 00237 

. 00328 

. 00427 

. 00560 

. 00702 

00861 

. 01037 

. 01232 

. 01445 

.01676 

. 01926 

. 02186 

. 02482 

02789 

.03116 

. 03462 

. 03828 

. 04214 

. 04621 

. 05049 

. 05495 

05970 

. 06463 

. 06978 

. 07517 

. 08077 

. 08664 

. 09274 

. 09908 

10555 

. 11248 

. 11955 

. 12686 

. 13442 

. 14218 

. 15029 

.15860 

16716 

. 17582 

. 18506 

. 19439 

. 20399 

. £1386 

. 22404 

. £3439 

£4505 

. £5598 

. £6729 

. £7859 

. £9026 

. 30217 

. 31431 

. 32645 

33928 

. 35209 

. 36511 

. 37859 

. 39176 

. 40538 

. 41919 

. 43317 

44754 

. 46166 

. 47615 

. 49079 

. 50577 

. 52053 

. 53561 

. 55081 

56615 

. 58243 

.59716 

.61283 

. 62859 

. 64516 

. 66038 

. 67675 

69355 

. 71036 

. 72667 

. 74406 

. 76094 

. 77783 

. 79284 

.81169 

82868 

. 84572 

. 86276 

. 87989 

. 89702 

. 91416 

. 93131 

. 94859 

96565 

. 98282 

1 . 00000 








THE NATURAL FREQUENCY OF VIBRATION (RAD/SEC) 
*****************■•#•*•*•****•#■**■**■#■**•****•■**•*■*•*•*■** 

SB . 06614 


THE NORMALIZED DISPLACEMENTS AT THE LUMPED MASSES 
************************************************* 


00000 

00001 

00005 

00032 

-. 00090 

-. 00135 

-. 00151 

-.00211 

00370 

00617 

00953 

-.01375 

-. 01884 

-.02428 

-.03153 

-.03912 

04752 

05671 

06669 

-. 07742 

-. 08889 

-. 10108 

-. 11355 

-. 12749 

14165 

15642 

17177 

-. 18766 

-. 20406 

-. 22093 

-. 23824 

-.25583 

27403 

29242 

-.31109 

-. 32999 

-. 34903 

-. 36827 

-. 38754 

-. 40682 

42574 

44517 

-. 46412 

-. 48283 

-.50124 

-. 51917 

-. 53688 

-. 55397 

57048 

58605 

-. 60143 

-. 61571 

-. 62907 

-. 64144 

-. 65277 

-. 66281 

67166 

67920 

-. 68542 

-. 69010 

-. 69335 

-. 69507 

-. 69520 

-.69374 

69051 

68561 

-. 67894 

-. 67031 

-. 66020 

-. 64806 

-. 63404 

-.61811 

59997 

58045 

-. 55869 

-. 53494 

-. 50887 

-. 48144 

-. 45167 

-.41989 

38610 

34836 

-. 31255 

-. 27283 

-. 23118 

-. 18563 

-. 14227 

-. 09425 

04430 

. 00634 

. 05600 

. 10952 

. 16198 

. 21499 

. 26256 

. 32365 

38036 

. 43859 

. 49804 

. 55884 

. 62045 

. 68279 

. 74569 

. 80940 

87253 

. 93623 

1 . 00000 








THE NATURAL FREQUENCY OF VIBRATION (RAD/SEC) 
******************************************** 

158. 90858 


THE NORMALIZED DISPLACEMENTS AT THE LUMPED MASSES 
**************•*•****■#•#*■***•**•*•*•**■***■#*■*•*■**■•*■*•#■**•*■*•** 


00000 

. 00001 

. 00010 

. 00069 

. 00190 

. 00286 

. 00318 

. 00443 

00770 

. 01270 

. 01938 

. 02768 

. 03752 

. 04790 

. 06153 

. 07556 

09082 

. 10722 

. 12469 

. 14310 

. 16237 

. 18237 

. 20237 

. 22414 

24567 

. 26746 

. 28940 

.31 136 

. 33320 

. 35478 

. 37598 

. 39652 

41666 

. 43585 

. 45409 

. 47121 

. 48705 

. 50154 

. 51448 

. 52579 

53520 

. 54303 

. 54878 

. 55249 

. 55408 

. 55348 

. 55061 

. 54544 

53791 

. 52818 

. 51564 

. 50088 

. 48368 

. 46406 

. 44193 

.41772 

39118 

. 36257 

. 33169 

. 29978 

. 26595 

. 23077 

. 19443 

. 15788 

11924 

. 08088 

. 04235 

.00319 

~. 03412 

07149 

-. 10790 

14303 

17707 

20825 

23770 

-. 26461 

-. 28891 

30946 

-. 32672 

-. 34010 

34927 

35402 

35369 

-. 34833 

-. 33755 

32021 

-. 29876 

-. 27061 

23888 

20453 

16881 

-. 12819 

— . 08630 

-.04199 

-. 00032 

. 05892 

12062 

. 18977 

. 26528 

. 34671 

. 43267 

. 52241 

. 61504 

.71040 

80585 

. 9027 1 

1 . 00000 








THE NUMBER OF NATURAL FREQUENCIES AT THE LISTED ROTOR SPEED (RPM) 
***************************************************************** 


3 40. 0 


THE NATURAL FREQUENCY OF VIBRATION (RAD/SEC) 
******************************************** 

11.52035 


THE NORMALIZED DISPLACEMENTS AT THE LUMPED MASSES 
************************************************* 


. 00000 

. 00000 

. 00001 

. 00006 

. 00065 

.00110 

. 00171 

. 00249 

. 00901 

. 01085 

. 01288 

. 01510 

. 02906 

. 03244 

. 03602 

. 03981 

. 06189 

. 06696 

. 07225 

. 07778 

. 10887 

. 11594 

. 12314 

. 13059 

. 17151 

. 18028 

. 18962 

. 19906 

. 25012 

.26110 

. 27246 

. 28380 

. 34458 

. 35738 

. 37038 

. 38384 

. 45253 

. 46658 

. 48099 

. 49554 

. 57035 

. 58650 

. 60111 

. 61664 

. 69660 

.71326 

. 72941 

. 74663 

. 83042 

. 84729 

. 86416 

. 88112 

. 96600 

. 98300 

1 . 00000 



00016 

. 00023 

. 00026 

. 00037 

00345 

. 00448 

. 00587 

. 00735 

01750 

. 02010 

. 02280 

. 02588 

04380 

. 04800 

. 05241 

. 05701 

08353 

. 08954 

. 09578 

. 10226 

13828 

. 14617 

. 15440 

. 16283 

20875 

. 21870 

. 22897 

. 23939 

29551 

. 30744 

. 31960 

. 33175 

39697 

. 41055 

. 42430 

. 43823 

51042 

. 52508 

. 54005 

. 55514 

63226 

. 64868 

. 66375 

. 67997 

76334 

. 78007 

. 79493 

. 81360 

89807 

. 91503 

. 93201 

. 94912 



THE NATURAL FREQUENCY OF VIBRATION <RAD/5EC> 
******************************************** 

69. 19415 


THE NORMALIZED DISPLACEMENTS AT THE LUMPED MASSES 
■************************************************* 


00000 

00001 

00005 

00033 

00090 

00136 

00151 

00212 

00371 

00619 

00955 

01378 

01887 

02431 

03157 

03916 

04756 

05675 

06671 

07743 

08889 

10105 

11350 

12739 

14152 

15625 

17154 

18737 

20370 

22049 

23772 

25522 

27332 

29161 

31016 

32893 

34784 

36694 

- . 38606 

40519 

42396 

44322 

46200 

48053 

49876 

51651 

53404 

55094 

56726 

58265 

59784 

61194 

62514 

— . 63734 

64851 

65839 

66709 

67450 

68060 

68517 

68833 

68997 

— . 69003 

68852 

68526 

68034 

67368 

66506 

65499 

64291 

62895 

6131 1 

59508 

57567 

55404 

53044 

50453 

47728 

44770 

41611 

38252 

34499 

30939 

26988 

22844 

18310 

13993 

09212 

04237 

. 00808 

. 05756 

. 11089 

. 16318 

. 21602 

. 26344 

. 32437 

38094 

. 43907 

. 49843 

. 55915 

. 62070 

. 68298 

. 74584 

. 80951 

87260 

. 93626 

1 . 00000 








THE NATURAL FREQUENCY OF VIBRATION (RAD/SEC) 
***********•*■*■**#■**#•*•*■•*****•*■**■#**•**■**•*■**#■**** 

159. 80793 


THE NORMALIZED DISPLACEMENTS AT THE LUMPED MASSES 
*##*#*###**#*#*#************#******###***#####**# 


00000 

. 00001 

. 00010 

. 00069 

. 00190 

. 002 Q 6 

. 00319 

. 00444 

00772 

.01272 

. 01941 

. 02771 

. 03756 

. 04794 

. 06158 

.07561 

09086 

. 10726 

. 12472 

. 14312 

. 16237 

. 18236 

. 20233 

. 22407 

£4556 

. £6732 

. 28921 

.31112 

. 33290 

. 35443 

. 37558 

. 39605 

41613 

. 43526 

. 45342 

. 47048 

. 48625 

. 50068 

. 51356 

. 52480 

53415 

. 54194 

. 54764 

. 55131 

. 55286 

. 55225 

. 54936 

. 54418 

53666 

. 52695 

. 51444 

. 49971 

. 48257 

. 46302 

. 44097 

.41684 

39040 

. 36190 

. 331 14 

. £9935 

. £6566 

. £3060 

. 19440 

. 15799 

11949 

. 08126 

. 04285 

. 00381 

-. 03339 

-. 07067 

-. 10698 

-. 14204 

17603 

-. £0718 

-. £3661 

-. 26352 

-. £8784 

-. 30843 

-. 32575 

-. 33921 

34848 

-. 35335 

-. 35314 

-. 34792 

-. 33728 

-. 32010 

29879 

-. £7079 

23920 

-. £0499 

-. 16940 

-. 12890 

-. 08713 

-. 04291 

-. 00132 

. 05785 

11953 

. 18869 

. £6425 

. 34576 

. 43181 

. 52166 

.61444 

. 70994 

80554 

. 90255 

1 . 00000 








4.6 EXAMPLE OF MODE SHAPE PLOTTING 


The following example illustrates the use of the MODPLT 
program to plot the mode shapes on a Tektronix 4010 or 4014 
display terminal. The program is executed by entering MODPLT. 
The program then queries the user to the computer interface 
connected to the terminal and the type of Tektronix terminal being 
used. The program then allows the user to generate a sketch of 
the mode shape; to change the scale factor used in sketching the 
mode shape; to make a complete finished drawing of the mode shape 
or to go on to the mode shape corresponding to the next natural 
frequency. For each mode shape, the normal procedure is to 
generate several sketches in order to determine a suitable scale 
factor and then generate a finished drawing which includes the 
title of the problem, the mode of vibration, the corresponding 
natural frequency, the number of lumped masses in the model, the 
rotor speed and the current date. 



A> MODPLT 

ENTER THE DATE (DD-MMM-YY) 

EB-DEC-Q5 

ENTER TYPE OF DISPLAY TERMINAL (4010 OR 4014) 

4014 

1 = DRAW SKETCH 

£ = CHANGE AMPLITUDE SCALE FACTOR, CURRENTLY = 

3 = DRAW FINISHED DRAWING (S) 

4 = CONSIDER NEXT NATURAL FREQUENCY AND MODE SHAPE 
1 




1 = DRAW SKETCH 

2 = CHANGE AMPLITUDE SCALE FACTOR, CURRENTLY = .50 

3 = DRAW FINISHED DRAWING <S) 

A = CONSIDER NEXT NATURAL FREQUENCY AND MODE SHAPE 
2 

ENTER NEW SCALE FACTOR TO BE APPLIED TO AMPLITUDE 
■ 70 

1 = DRAW SKETCH 

on 2 = CHANGE AMPLITUDE SCALE FACTOR, CURRENTLY = .70 

w 3 = DRAW FINISHED DRAWING (S) 

A = CONSIDER NEXT NATURAL FREQUENCY AND MODE SHAPE 
1 




1 = DRAW SKETCH 

£ ■ CHANGE AMPLITUDE SCALE FACTOR, CURRENTLY = .70 

3 = DRAW FINISHED DRAWING <S) 

4 = CONSIDER NEXT NATURAL FREQUENCY AND MODE SHAPE 


a\ 

Ui 



FLRPWISE VIBRRTION OF THE BLRDE 



MODE 

NAT. FREQ. 
<HZ.) 

NO. OF RASSES 

ROTOR 

RPO 

PATE 

1 

1.646 

107 

.0 

28-CEC-85 




1 = DRAW SKETCH 

£ = CHANGE AMPLITUDE SCALE FACTOR, CURRENTLY = .70 

3 = DRAW FINISHED DRAWING <S> 

4 = CONSIDER NEXT NATURAL FREQUENCY AND MODE SHAPE 
4 

T = DRAW SKETCH 

£ = CHANGE AMPLITUDE SCALE FACTOR, CURRENTLY = .70 

^ 3 = DRAW FINISHED DRAWING (S> 
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NODE NAT • FREQ. NO. OF HASSES ROTOR DATE 

(HZ.) RPH 

2 It. 773 107 .• 28-DEC-8S 




FLRPWISE VIBRRTION OF THE BLRDE 
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BODE NAT. FREO. NO. OF HASSES ROTOR DATE 

(HZ. > RPN 

3 25.243 107 .0 28-DEC-8S 




FLRPWISE VIBRRTION OF THE BLRDE 
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BODE NAT. FREO. NO. OF BASSES ROTOR PATE 

(HZ. > RPfj 

1 1.69S 1#7 SO.* 88-PEC-85 




FLAPWISE VIBRATION OF THE BLADE 
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IT. FREQ. NO. OF BASSES ROTOR DATE 

<HZ. ) RPH 

IQ. 833 1Q7 20.0 28-0EC-8S 
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NODE NAT. FREO. HO. OF HASSES ROTOR DATE 
CH 2 . > RPH 

3 as.aai 107 aa.e 28-dec-ss 
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BODE MAT. FREQ. NO. OF BASSES ROTOR DATE 

(HZ. ) RPN 

1 1.834 107 49. 0, 28-DEC-85 
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MODE MOT. FREO. MO. OF MASSES ROTOR DATE 

(HZ. ) RPM 

2 11.013 107 40.0 28-DEC-8S 
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NAT. FREQ. NO. OF MASSES ROTOR DATE 

(HZ.) RPH 

as . 434 107 46. 6 88-DEC-85 




5. ANALYTICAL AND EXPERIMENTAL RESULTS 


The sections which follow document the analytical and 
experimental results for the blade and counterweight assemblies. 
The analytical results were obtained by executing the computer 
codes which are described in Sections 3 and 4. The experimental 
results were obtained using real time spectral analysis 
techniques . 


5.1 ANALYTICAL RESULTS FOR THE MOD-O BLADE ASSEMBLY 


The BLDFLAP computer code was used to generate a 573 lumped 
mass model of the blade assembly for flapwise vibration. This 
assembly included a rigid hub having a radius of 22.0 inches, the 
spool piece (CF 764254), the "Shorty 40" blade, the transition 
section and a six foot pitchable tip section. The number of 
lumped masses for each section was : 


Spool Piece Inboard Flange - 2 
Spool Piece Body - 14 
Spool Piece Outboard Flange - 2 
Blade Section - 440 
Transition and Tip Sections - 115 


This lumped mass model was used as input to the BEAM code and 
results were obtained for rotor speeds ranging from 0 rpm (static 
rotor condition) to 90 rpm at intervals of 10 rpm. Additional 
input to the BEAM program was: 


Lower Limit of Searching Frequency Range 
Upper Limit of Searching Frequency Range 
Searching Frequency Interval 
Number of Digits of Accuracy 


0.0 rad/sec 
350.0 rad/sec 
10.0 rad/sec 
3 


Natural frequencies and mode shapes were obtained for the first 
four modes of vibration. The natural frequencies for each of the 
rotor speeds are listed in Table 5.1. 


It can be seen that, as predicted, the natural frequencies 
increase with increasing rotor speed. This effect, although not 
extremely significant, is less pronounced at the higher modes of 
vibration . 
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TABLE 5.1 


THE FIRST FOUR NATURAL FREQUENCIES FOR FLAPWISE VIBRATION OF 
THE BLADE ASSEMBLY (CPS) AT VARIOUS ROTOR SPEEDS 


0 RPM 

10 RPM 

20 RPM 

30 RPM 

40 RPM 

1.64 

1.65 

1.69 

1.75 

1.83 

10.70 

10.72 

10.76 

10.84 

10.94 

25.02 

25.04 

25.07 

25.13 

25.21 

46.13 

46.14 

46.18 

46.24 

46.33 


50 RPM 

60 RPM 

70 RPM 

80 RPM 

90 RPM 

1.92 

2.03 

2.16 

2.29 

2.43 

11 .07 

11.23 

11.42 

11 .63 

11.86 

25.32 

25.45 

25.60 

25.77 

25.97 

46.45 

46.59 

46.76 

46.95 

47.16 
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5.2 ANALYTICAL RESULTS FOR THE COUNTERWEIGHT ASSEMBLY 


The CWTFLAP computer code was used to generate a 403 lumped 
mass model of the counterweight assembly for flapwise vibration. 
This assembly included a rigid hub having a radius of 22.0 inches, 
the spool piece (CF 764254), the tapered steel spar extension (CF 
760549), and the counterweight (CF 764554). The number of lumped 
masses for each section was : 


Spool Piece Inboard Flange - 6 
Spool Piece Body - 42 
Spool Piece Outboard Flange - 6 
Spar Base Flange - 6 
Spar Inboard Transition Portion - 15 
Spar Tapered Portion - 309 
Spar Outboard Transition Portion - 12 
Spar Tip Flange - 6 
Lumped Rigid Counterweight Mass - 1 


The lumped weight of the rigid counterweight mass located at the 
free end of the tapered spar extension was 5294 pounds . This 
lumped mass model was used as input to the BEAM code and results 
were obtained for rotor speeds ranging from 0 rpm {static rotor 
condition) to 90 rpm at intervals of 10 rpm. Additional input to 
the BEAM program was : 


Lower Limit of Searching Frequency Range 
Upper Limit of Searching Frequency Range 
Searching Frequency Interval 
Number of Digits of Accuracy 


0.0 rad/sec 
7000.0 rad/sec 
100.0 rad/sec 
3 


Natural frequencies and mode shapes were obtained for the first 
four modes of vibration. The natural frequencies for each of the 
rotor speeds are listed in Table 5.2. 


Although the natural frequencies of the counterweight assembly 
increase with increasing rotor speed, the effect is even less 
pronounced than that of the blade assembly. An identical lumped 
mass model of the counterweight assembly was used to perform a 
chordwise vibration analysis for the static rotor condition. The 
natural frequencies for the first four modes of vibration are 
listed in Table 5.3. It should be noted that for static rotor 
conditions , a chordwise and flapwise vibration analysis of the 
counterweight assembly yielded almost identical results as can be 
seen in Tables 5.2 and 5.3. 
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TABLE 5.2 


THE FIRST FOUR NATURAL FREQUENCIES FOR FLAPWISE VIBRATION OF 
THE COUNTERWEIGHT ASSEMBLY (CPS) AT VARIOUS ROTOR SPEEDS 


0 RPM 

10 RPM 

20 RPM 

30 RPM 

40 RPM 

8.21 

8.21 

8.22 

8.23 

8.25 

171.11 

171.12 

171.14 

171.17 

171.21 

526.11 

526.12 

526.14 

526.17 

526.22 

1034.82 

1034.83 

1034.85 

1034.88 

1034.93 


50 RPM 

60 RPM 

70 RPM 

80 RPM 

90 RPM 

8.27 

8.30 

8.33 

8.37 

8.41 

171.27 

171.33 

171.41 

171.50 

171.61 

526.28 

526.35 

526.44 

526.54 

526.65 

1034.99 

1035.06 

1035.15 

1035.25 

1035.37 
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TABLE 5.3 


THE FIRST FOUR NATURAL FREQUENCIES FOR CHORDWISE VIBRATION OF 
THE COUNTERWEIGHT ASSEMBLY (CPS) AT A ROTOR SPEED OF ZERO 


8.21 

171.12 

526.16 

1035.00 
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5.3 EXPERIMENTAL RESULTS FOR THE BLADE ASSEMBLY 


Frequency domain spectrum were supplied by NASA personnel for 
a modal analysis conducted at the Plum Brook facility on the blade 
assembly. The tip sections on the two-bladed rotor assembly were 
removed and a modal analysis was performed on the resulting 
assembly in a teetered hub configuration for both chordwise and 
flapwise vibration. For flapwise vibration, the response spectrum 
showed peaks at frequencies of 2.92 cps , 9.0 cps , 13.3 cps , 16.8 
cps , 25.3 cps, 30.5 cps, and 39.0 cps. The analytical values for 
the first four modes of vibration of a model including the tip 
section are given in Table 5.1 as 1.64 cps, 10.70 cps, 25.02 cps 
and 46.13 cps. The blade model generation program, BLDFLAP , was 
modified to remove the tip section. As a result of removing the 
tip section the natural frequencies were found to increase. The 
analytical values for the first four modes of vibration of the 
blade assembly without the tip section were found to be 1.94 cps, 
13.14 cps, 33.93 cps and 64.83 cps. The additional peaks found 
experimentally may be due to the rigid body motion of the assembly 
and the effects of the rubber stops or other nonlinear effects in 
the system. 


5.4 EXPERIMENTAL RESULTS FOR THE COUNTERWEIGHT ASSEMBLY 


Two trips were made to the Lewis Research Center to acquire 
experimental modal analysis data for the counterweight assembly. 
This data was obtained using both impact excitation techniques 
which gave the structure an initial velocity without an initial 
displacement and displacement techniques which resulted in giving 
the structure an initial displacement without an initial velocity. 
The response of the structure was obtained using both an 
accelerometer and a microphone. Output from these sensors was 
processed using a single channel Spectral Dynamics Real Time 
Analyzer and the results were recorded using a Hewlett Packard X-Y 
Plotter. All equipment was provided by the Mechanical Engineering 
Department of The University of Toledo. The counterweight 
assembly including the spool piece was mounted on a test fixture 
supplied by NASA. Data was taken for flapwise vibration of the 
counterweight assembly for static rotor conditions. 


Typical analysis output (frequency domain spectra) are shown 
in Figure 5.1 for an analysis range of 20 cps and in Figure 5.2 
for an analysis range of 2000 cps . The vertical axis represents 
the amplitude of the signal in terms of a logarithmic scale and 
the horizontal axis represents the frequency of the spectra in 
terms of a linear scale. The first peak shown in Figure 5.1 
occurs at 1.44 cps and represents internal noise generated within 
the real time analyzer. The second peak occurs at 6.36 cps and is 
believed to be the fundamental frequency of the counterweight 
assembly as mounted on the test fixture. A model of the 
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Figure 5.2 Frequency Domain Spectra for Counterweight Assembly 
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counterweight assembly which did not include the spool piece {the 
spool piece was originally considered to be rigid with respect to 
the rest of the assembly) yielded a fundamental frequency of 9.09 
cps for flapwise vibration. The counterweight model generation 
program, CWTFLAP , was then modified to include the effects of the 
spool piece. The model which included the spool piece yielded an 
analytical value for the fundamental frequency of 8.21 cps for 
flapwise vibration of the counterweight assembly. The discrepancy 
between the experimental value of 6.36 cps and the analytical 
value of 8.21 cps for the flapwise fundamental frequency of the 
counterweight assembly may very well be due to the lack of 
rigidity in the test fixture which supported the assembly. In 
order to determine if this was indeed the case, one would need to 
perform a complete modal analysis of the counterweight and test 
fixture assembly using a two or four channel modal analysis system 
which is capable of generating animated mode shapes as well as 
performing a frequency analysis. Such portable equipment was not 
readily available at the time this study was conducted. 


The model of the counterweight which included the spool piece 
predicted frequencies of 171.11 cps, 526.11 cps and 1034.82 cps 
for the second, third and fourth modes of vibration. As seen in 
Figure 5.2 there were many peaks that were experimentally acquired 
in this range. In fact, there were so many peaks that it was 
impossible to experimentally determine the remaining natural 
frequencies of the counterweight assembly. It is believed that 
these peaks resulted from the ringing effects and the pipe organ 
effects of the thin shelled steel spar extension. As predicted by 
the analytical mode shapes and noted in the field, the free end of 
the assembly where the counterweight is located became a nodal 
point (a point of zero displacement) for vibration at natural 
frequencies above the fundamental mode. Thus it was necessary to 
attach the accelerometer directly to the steel spar extension when 
trying to determine the natural frequencies for the higher modes 
of vibration. It is believed that this resulted in picking up not 
only the vibration of the assembly as a cantilevered beam but the 
radial vibration of the steel spar extension as a thin shell of 
revolution as well as the vibration of the column of air trapped 
inside the assembly. Also, the question of the interaction 
between the counterweight assembly and the non-rigid test fixture 
remains unanswered. In order to determine the natural frequencies 
for flapwise vibration of the counterweight assembly from Figure 
5.2, one would need to use a multi-channel modal analysis system. 
Such a system would not only determine the frequency response 
functions but would determine the phase lag between the system 
response and the excitation. This additional information would 
allow the determination of the animated mode shapes for each of 
the peaks shown in Figure 5.2 and would thus yield the natural 
frequencies for flapwise vibration of the assembly. 
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DISCUSSION OF RESULTS AND CONCLUSIONS 


6 . 


The sections which follow discuss the results and conclusions 
of this study. The results and conclusions of the analysis of the 
blade assembly and the counterweight assembly are discussed in 
Sections 6.1 and 6.2 respectively. Section 6.3 compares the 
application of a personal computer to perform the vibration 
analysis with the application of a minicomputer. Recommendations 
for further study are given in Section 6.4. 


6.1 DISCUSSION OF THE MOD-O BLADE ASSEMBLY 


It should be emphasized that the discussion which follows is 
based on the uncoupled vibration of the blade assembly treated as 
a cantilevered beam. The results presented in Section 5.1 
indicate that the natural frequencies for the flapwise vibration 
of the Mod-0 blade assembly increase with increasing rotor speeds. 
However, the percentage increase was not as significant at the 
higher natural frequencies. The natural frequencies for the first 
four modes of vibration were significantly higher than the rotor 
speed which was permitted to vary between 0 and 90 rpm. However, 
considering the rotor is excited twice per revolution, the 
fundamental natural frequency can be excited at a rotor speed 
between 50 and 60 rpm. Since the flexural rigidity of the blade 
assembly is significantly stiffer in the chordwise direction, the 
natural frequencies are above the operating range and a chordwise 
analysis was not performed. 


6.2 DISCUSSION OF THE COUNTERWEIGHT ASSEMBLY 


As in the case of the blade assembly, it should be emphasized 
that the discussion which follows is based on the uncoupled 
lateral vibration of the counterweight assembly treated as a 
cantilevered beam. The results presented in Section 5.2 indicate 
that the natural frequencies for the flapwise vibration of the 
counterweight assembly increase with increasing rotor speeds 
although this effect is even less pronounced than that with the 
blade assembly. The percentage increase was not as significant at 
the higher natural frequencies . Natural frequencies for the first 
four modes of vibration were significantly higher than the rotor 
speed which was permitted to vary between 0 and 90 rpm. Even, 
considering the rotor is excited twice per revolution, the 
fundamental natural frequency remains well above the forcing 
frequency. Since the flexural rigidity of the counterweight 
assembly is essentially the same in both the flapwise and 
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chordwise directions , the natural frequencies for vibration in the 
chordwise direction are almost identical to those for flapwise 
vibration . 


6.3 COMPARISON OF PC AND MINICOMPUTER COMPUTATIONS 


One of the primary goals of this research effort was to permit 
users to employ a personal computer for the vibratory analysis of 
wind turbines. Traditionally this type of analysis was restricted 
to mainframe and minicomputer environments. The analysis codes 
used in this research effort were originally developed on a DEC 
PDP 11/34 minicomputer. This code was then downloaded and 
modified to run on an ITT XTRA , IBM compatible, personal computer. 
Studies were conducted to evaluate the efficiency of running the 
wind energy programs on a personal computer as compared to the 
minicomputer. These studies were conducted using a 139 mass model 
of the counterweight assembly in flapwise vibration at static 
rotor conditions. In all cases, the first four modes of vibration 
were determined. The execution times for various cases follow. 


System and Options 


Time ( secs ) 


PDP 11/34 with floating point hardware 715 

ITT PC with default math library but without 

8087 math co-processor 6031 

ITT PC with 8087 math co-processor and 

default math library 955 

ITT PC with 8087 math co-processor and 

8087 math library 815 

ITT PC with 8087 math co-processor, 8087 math 
library, and SNOFLOATCALLS and $STORAGE:2 

metacommands 644 


One can note that, with the proper combination of hardware and 
software options, the efficiency of using a PC exceeds that of a 
minicomputer. It is also noteworthy that, since the quoted 
studies were conducted, an additional hardware option has been 
released which further increases the speed of the ITT PC. Thus, 
the application of personal computers to this type of vibration 
analysis is very feasible. 
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6.4 RECOMMENDATIONS FOR FURTHER STUDY 


There are at least five additional areas of research which 
should be investigated in order to guarantee successful operation 
of a wind turbine rotor assembly. The first research area 
involves the development of a model generation program for the 
chordwise vibration of the blade assembly. Currently, models may 
only be generated for flapwise vibration of the blade assembly. 
The second research area involves the effect of rotation on the 
natural frequencies and mode shapes. The BEAM analysis code 
presently evaluates the effects of rotation for the flapwise 
vibration of the uncoupled blade and counterweight assemblies. 
This code should be modified to include the effects of rotation on 
chordwise vibration. The third research area involves the 
coupling of the blade and counterweight assemblies into a single 
rotor configuration. Such a configuration would require that the 
BEAM program be modified to reflect the altered boundary 
conditions . The changes are required because the previously 
assumed cantilevered beam model will have to be supplemented with 
a centrally supported teetered hub model. The fourth area of 
research involves the analysis of the tower shadow effects on the 
forcing frequency applied to the rotor assembly. Applying Fourier 
Series Analysis techniques to the tower shadow effect may reveal 
forcing frequencies which differ from the one or two per 
revolution that were assumed in the current research effort. 
Lastly, in an effort to obtain more definitive experimental 
results, a complete modal analysis of the blade and counterweight 
assemblies should be performed using a multichannel modal analysis 
system capable of generating animated mode shapes as well as 
performing a frequency analysis. Because of the various 
nonlinearities in the entire system, the mode shapes are necessary 
to properly identify the natural frequencies of vibration. 
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APPENDIX 


This appendix contains the listings of all computer codes 
described in Section 3. The order of the listings and their 
functions are: 

BLDFLAP * blade assembly model generation code 
CWTFLAP * counterweight assembly model generation code 
CWTCHRD *• counterweight assembly model generation code 
BEAM * main analysis code 

DET - subroutine to evaluate determinant 
FILLA - subroutine to fill the A matrix 
FORCE “ subroutine to calculate centrifugal forces 
INPUT •* subroutine to perform all input functions 
MODES + subroutine to calculate mode shapes 
MULT ** subroutine to perform matrix multiplication 
NFREQ - subroutine to determine natural frequencies 
OUTPUT - subroutine to perform all output functions 
MODPLT * mode shape plotting code 

DRAW * subroutine to draw mode shapes 

MINPUT - subroutine to perform all input functions. 
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^STORAGE : £ 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


A. 1 BLDFLAP •* i blade assembly model generation code 

THIS PROGRAM GENERATES THE INPUT DATA FOR THE SPOOL 
PIECE (CF 764254) , THE "SHORTY 40" BLADE, TRANSITION 
SECTION AND THE S’ PITCHABLE TIP FOR FLAPWISE VIBRATION 
(VIBRATION PERPENDICULAR TO THE PLANE OF ROTATION). THE 
PROGRAM CONSIDERS THE ASSEMBLY TO BE BROKEN UP INTO 4 
SECTIONS AS FOLLOWS: 


SECTION 1 - THE FLANGE AT THE BASE OF THE SPOOL PIECE 

OF LENGTH 2.25" WEIGHING APPROX. 10B. 74 LBS. 

SECTION 2 - THE WEBBED MAIN BODY OF THE SPOOL PIECE 

OF LENGTH 13.75" WEIGHING APPROX. 207.96 LBS. 

SECTION 3 - THE FLANGE AT THE SPAR END OF THE SPOOL PIECE 
OF LENGTH 1.75" WEIGHING APPROX. B2. 78 LBS. 

SECTION 4 - THE BLADE SECTION OF LENGTH 441.25" 

WEIGHING APPROX. 2250.00 LBS. 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SECTION 5 - THE TRANSITION AND TIP SECTIONS OF LENGTH 
115.56" WEIGHING APPROX. 960.00 LBS. 


El (N> THE PRODUCT OF MODULUS OF ELASTICITY AND AREA MOMENT 
OF INERTIA AT MASS N (LBS* INCH* INCH) 

Eli ( J) THE INPUT VALUES OF El USED TO DEFINE THE RIGIDITY OF 
THE BLADE, TRANSITION AND TIP SECTIONS. THESE VALUES 
ARE ENTERED USING A DATA STATEMENT IN THE PROGRAM. 
(LBS* INCH* INCH) 

EIO THE PRODUCT OF MODULUS OF ELASTICITY AND AREA MOMENT 

OF INERTIA AT ZERO OR THE FIXED END <LBS*INCH*INCH) 

INEI THE INPUT NUMBER OF El VALUES USED TO DEFINE THE 
RIGIDITY OF THE BLADE, TRANSITION SECTION AND TIP 
SECTION 

INW THE INPUT NUMBER OF VALUES OF THE WEIGHT BETWEEN 

STATIONS USED TO DEFINE THE WEIGHT DISTRIBUTION OF 
THE BLADE, TRANSITION SECTION AND TIP SECTION 

LEN(N) THE LENGTH OF THE SECTION TO THE LEFT OF MASS N 
(INCHES) 

NB THE MASS NUMBER AT THE BEGINNING OF A SECTION 

NF THE MASS NUMBER AT THE END OF A SECTION 

ND THE NUMBER OF DIVISIONS IN A PARTICULAR SECTION 
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n n n n n 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


NS THE NUMBER OF THE SECTION UNDER CONSIDERATION 

NTOT THE TOTAL NUMBER OF MASSES 

RRHUB THE RADIUS OF THE RIGID HUB (INCHES) 

TITLE THE PROBLEM TITLE STATEMENT 

WKL) THE INPUT VALUES OF THE WEIGHT BETWEEN STATIONS 
USED TO DEFINE THE WEIGHT DISTRIBUTION OF THE 
BLADE, TRANSITION AND TIP SECTIONS. THESE VALUES 
ARE ENTERED USING A DATA STATEMENT IN THE 
PROGRAM. (LBS. ) 

WPI (M) THE WEIGHT OF THE BuADE PER UNIT LENGTH CALCULATED 
FROM THE INPUT VALUES OF WKL) (LBS. /INCH) 

WT (N) THE WEIGHT OF MASS N (LBS) 

X (N) THE AXIAL LOCATION OF MASS N (INCH) 

XEII(J) THE AXIAL LOCATION WHERE THE INPUT VALUES OF El 
ARE SPECIFIED IN ORDER TO DEFINE THE RIGIDITY 
OF THE BLADE, TRANSITION AND TIP SECTIONS. THESE 
VALUES ARE ENTERED USING A DATA STATEMENT IN THE 
PROGRAM. (INCH) 

XWI(I) THE AXIAL LOCATIONS WHERE THE INPUT VALUES OF THE 
WEIGHT BETWEEN STATIONS ARE SPECIFIED IN ORDER TO 
DEFINE THE WEIGHT DISTRIBUTION OF THE BLADE, 

TRANSITION AND TIP SECTIONS. THESE VALUES ARE ENTERED 
USING A DATA STATEMENT IN THE PROGRAM. (INCH) 


IMPLICIT REALMS (A-H, O-Z) 

INTEGER R, RM1 
REAL*8 INTR, LEN 

DIMENSION El (700) , LEN (700) , TITLE (15) , WT(700) , X (700) , 
XEII (21) , XWI (25) , XWPI (30) , Eli (21) , WI (25) , WPI (30) 

INPUT BLADE, TRANSITION SECTION AND TIP SECTION 
PROPERTIES FOR THE WEIGHT DISTRIBUTION AND 
FLEXURAL RIGIDITY 


DATA INW/25/ 

DATA WI/1 15. 1 3 DO, 290. 80D0, 154. 81D0, 154. 36DO, 

1 121. 1 1 DO, 120. 02DO, 113. 64D0, 105. 34D0, 92. 21D0, 

2 82. 49D0, 70. 16D0, 75. 44D0, 84. 20DO, 69. 62D0, 85. 07D0, 

3 99. 1 2D0, 1 18. 36D0, 198. 80D0, 104. 32D0, 326. ODO, 

4 321 . ODO, 98. ODO, 91 . ODO, 84. ODO, O. ODO/ 

DATA XWI/39. 75D0, 50. ODO, 75. ODO, 100. ODO, 125. ODO, 150. ODO, 

1 175. ODO, 200. ODO, 225. ODO, 250. ODO, 275. ODO, 300. ODO, 325. ODO, 

2 350. ODO, 375. ODO, 400. ODO, 425. ODO, 450. ODO, 475. ODO, 481 . ODO, 

3 500. ODO, 525. ODO, 550. ODO, 575. ODO, 596. 56D0/ 
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DATA INEI/21/ 

DATA El 1/285. D8, 208. DS, 115. D8, 65. D8, 51. D8, 35. D8, 

1 93. 96D8, 84. 86D8, 82. 57D8, 73. 47D8, 64. 37D8, 19. 52D8, 19. 52D8, 

2 15. 32D8, 14. 63D8, 13. 95D8, 13. 26D8, 12. 58D8, 1 1. 89D8, 

3 11. 53D8, 0. 000/ 

DATA XEI I /39. 75D0, 100. DO, 200. DO, 300. DO, 400. DO, 480. 9D0, 

1 481. DO, 491. DO, 501. DO, 511. DO, 521. DO, 521. 001D0, 527. 1S7D0, 

2 527. 188D0, 537. 187D0, 547. 187D0, 557. 187D0, 567. 187D0, 577. 187D0, 

3 582. 437D0, 596. 56D0/ 

C 

C INTERPOLATION FUNCTION 

C 

INTR (XX, XXL, XXR, FFL, FFR) = (FFR-FFL) * ( XX-XXL) / < XXR-XXL) +FFL 
C 

C ASSIGN LOGICAL UNIT 1 TO ’VIBIN.DAT’ 

C 

0PEN(1, FILE =’ VIBIN. DAT’ , STATUS =’ NEW’ ) 

C 

C READ DATA FROM TERMINAL 

C 

WRITE!*, 5) 

5 FORMAT!//, 2X, ’ THIS PROGRAM GENERATES THE DATA FOR THE ’, 

1 ’ FLAPWISE ’ , /, 2X, ’ VIBRATION OF THE BLADE ASSEMBLY ’, 

2 ///, £X , ’ ENTER THE RADIUS OF THE RIGID HUB ! INCHES) : ’ , <J>) 

READ!*, *) RRHUB 

WRITE!*, 10) 

10 FORMAT !//£X, ’ SECTION 1-SPOOL PIECE FLANGE L=2. 25" W=108.74 LB’, 

1 / /£X, ’ SECTION 2-SPOOL PIECE BODY L=13. 75" W=207. 96 LB’, 

2 / /£X, ’ SECTION 3-SPOOL PIECE FLANGE L=1.75" W=82. 78 LB’, 

4 //EX,’ SECTION 4-BLADE SECTION L=441.25" W=2£50. 00 LB’ 

5 //2X, ’ SECTION 5-TRANSITION & TIP SECTION L=115.56" W=960. OO LB’ X 

WRITE!*, £0) 

£0 FORMAT!/, 4X, ’ THE MAXIMUM TOTAL NUMBER OF DIVISIONS YOU ’ , 

1 / , 4X, ’ ARE ALLOWED IN ALL 5 SECTIONS IS 699. ’ ) 

E= 30000000. O 
PI=3. 141592654 
NS= 1 
NF=0 

30 GO TO !40, 60, 80, 100, 310, 400) NS 

C 

C SECTION 1 PARAMETERS 

C ******************** 

c 

40 WRITE!*, 50) 

50 FORMAT!//, ’ ENTER THE NUMBER OF DIVISIONS FOR SECTION 1 : ’ , $) 

READ !*, 104) ND 
WC=48. 3268 
XL=0. 

XR=2. 25 
D0=£ 1 . 63 
DI=15. 825 
C 

C CALCULATIONS FOR SECTIONS 1 AND 3 

C ******♦***#*#********#**####***** 

C 

55 DN=ND 
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DX= ( XR— XL) /DN 
XR=XL+DX 
NB=NF+1 
NF=NF+ND 

S£I=(E*PI*(D0**4-DI**4) ) /64. O 
DO 57 N=NB, NF 
WT (N) =WC* (XR— XL) 

X (N>=. 50* ( XR+XL ) 

El <N)=SEI 
XL=XR 

* 57 X R=XL+DX 

NS=NS+1 
GO TO 30 

1 C 

C SECTION £ PARAMETERS 

C ******************** 

c 

60 WRITE (*,70) 

70 FORMAT (//, ’ ENTER THE NUMBER OF DIVISIONS FOR SECTION £ : ’ , *> 

READ (*, 104) ND 
WC=15. 1£4£ 

XR=16. 0 

C 

C CALCULATIONS FOR SECTION £ 

C ************************** 

c 

5 DN=ND 

DX= <XR— XL) /DN 
XR=XL+DX 
NB=NF+1 
NF=NF+ND 
SEI=E*£106. £68 
DO 75 N=NB, NF 
WT(N)=WC*DX 
X (N) =. 50* (XR+XL) 

El <N)=SEI 
XL=X R 

75 XR=XL+DX 

NS=NS+1 
GO TO 30 

C 

C SECTION 3 PARAMETERS 

C ******************** 

c 

80 WRITE (*,90) 

90 FORMAT (//, ’ ENTER THE NUMBER OF DIVISIONS FOR SECTION 3 : ’ , *) 

READ (*, 104) ND 
WC=47. 30£1 
XR=17. 75 

* D0=£1 . 63 
DI=15. 97 
GO TO 55 

» C 

C SECTION 4 PARAMETERS 

C ******************** 
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c 

100 WRITE (*,1*02) 

102 FORMAT (//, ’ ENTER THE NUMBER OF DIVISIONS FOR SECTION 4 : ’ ,*) 

READ ( *, 104) ND 
104 FORMAT (15) 

XR=459. 00 
XL=17. 75 
C 

C INITIALIZE LEFT INDEX L. L WILL POINT TO THE WPI AND XWI VALUE 

C WHICH IS TO THE LEFT OF OR AT XL. 

C 

L=1 

C 

C CALCULATIONS FOR SECTIONS 4 AND 5 

C ********************************* 

C 

C CALCULATE THE X COORDINATES AS REFERENCED FROM RIBID HUB 

C 

DO 109 1=1, INW 

109 XWI (I)=XWI (I)-RRHUB 
DO 110 1=1, INEI 

110 XEII (I)=XEII (I)-RRHUB 
C 

C CALCULATE THE WT/IN FOR EACH GIVEN SECTION 

C 

INWM1=INW— 1 

DO 115 1=1, INWM1 

WPI <I)=WI (I) / (XWI (I + D-XWI (I) ) 

115 CONTINUE 
C 

C CALCULATE LUMPED WEIGHTS, AXIAL CG LOCATIONS AND El VALUES 

C 

116 DX=(XR-XL) /FLOAT (ND) 

C 

C CONSIDER ND MASSES 

C 

NB=NF+1 
NF=NF+ND 
DO 300 N=NB, NF 
XR=XL+DX 
C 

C FIND RIGHT INDEX R. R WILL POINT TO THE WPI AND XWI VALUE 

C WHICH CONTI ANS XR. 

C 

DO 120 R=L, INW 

IF(XR.GE. XWI (R). AND. XR.LE. XWI (R+l) ) GO TO 130 
120 CONTINUE 

R=INW 
C 

C CHECK IF DX IS TOTALLY WITHIN A WPI VALUE (ICASE=1>, IF DX 

C SPANS ACROSS TWO WPI VALUES <ICASE=2), OR IF DX SPANS 

C ACROSS 3 OR MORE VALUES (ICASE=3) 

C 

130 ICASE=R— L+l 

IF ( ICASE. GT. 3) ICASE=3 
GO TO (140, 180, ISO) , ICASE 
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c 

C CASE 1 - NEW MOSS ENTIRELY BETWEEN WPI VALUES 

C 

140 WT (N) =WPI (L) *DX 

X (N) = (XL+XR) /£. 

GO TO £40 
C 

C CASE £ AND 3 - NEW MASS SPANS TWO OR MORE WPI VALUES 

C 

180 RM1=R— 1 

LP1=L+1 

WTSUM=WPI (L) * (XWI (LPl)-XL) 

IF ( ICASE. EQ. £) GO TO 195 
DO 190 I=LP1 , RM1 

WTSUM=WTSUM+WPI (I)*(XWI <I+1)-XWI (I) ) 

190 CONTINUE 

C 

C ADD LAST PIECE 

C 

195 WTSUM=WTSUM+WP I (R)*(XR-XWI (R) ) 

C 

C CALCULATE WT*XCG FOR LEFT PARTIAL PORTION 

C 

XCGSUM=WPI (L) * < XWI (LP1 ) -XL) * (XWI <LP1 ) +XL) /£. 

C 

C ADD WT*XCG FOR RIGHT PARTIAL PORTION TO SUM 

C 

XCGSUM=XCGSUM+WPI <R)*<XR-XWI <R) )*(XR+XWI <R> ) /£. 

C 

C AVOID ADDING EXTRA TERMS TO SUM FOR CASE £ 

C 

IF ( ICASE. EQ.£) GO TO £10 
DO £00 I=LP1, RM1 
IP1=I+1 

XCGSUM=XCGSUM+WPI (I)*(XWI (I PI) -XWI (I) )*(XWI (IPD+XWI ( I ) ) /£. 
£00 CONTINUE 

£10 WT(N)=WTSUM 

X <N) =XCGSUM/WTSUM 
C 

C FIND El AT CALCULATED CG AXIAL LOCATION 

C 

240 INEIM1=INEI-1 

DO £50 LEI=1, INEIM1 

IF(X (N) . GE. XEII (LEI) . AND. X (N) . LT. XEII (LEI + 1) ) GO TO £60 
£50 CONTINUE 

LEI=NEIM1 

£60 El (N) =INTR ( X (N) , XEII (LEI ) , XEI I (LEI + 1 ) , El I (LEI ) , El I (LEI + 1 ) ) 

C 

C ADJUST VALUES FOR NEXT MASS 

C 

XL=XR 

DO £80 I=L, INW 

IF(XL.GE. XWI (I). AND. XL.LT. XWI (1+1) ) GO TO £90 
£80 CONTINUE 

L=INW 
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GO TO 300 
£90 L= I 

300 CONTINUE 

NS=NS+1 
GO TO 30 
C 

C SECTION 5 PARAMETERS 

C ******************** 

c 

310 WRITE (*,320) 

3£0 FORMAT (//, ’ ENTER THE NUMBER OF DIVISIONS FOR SECTION 5 : ’,*> 

READ (*, 330) ND 
330 FORMAT < 15) 

XR=574. 56 
XL=459. 00 
C 

C INITIALIZE LEFT INDEX L. L WILL POINT TO THE WPI AND XWI VALUE 

C WHICH IS TO THE LEFT OF OR AT XL. 

C 

L=20 

GO TO 116 
C 

C ENTER TITLE STATEMENT 

C 

400 WRITE (*, 410) 

410 FORMAT ( ///, ’ ENTER A TITLE STATEMENT (60 CHAR. MAX.) : ’,*) 

READ (*, 420) TITLE 
4£0 FORMAT (15A4) 

C 

C CALCULATE El AT THE BASE OR FIXED END OF THE SPOOL 

C PIECE (X=0. O) 

C 

D0=£1. 63 
DI=15. 8£5 

EI0=(E*PI*(D0**4-DI**4> > /64. O 
C 

C OUTPUT RESULTS TO VIBIN DATA FILE 

C ********************************* 

C 

700 WRITE <1, 710) TITLE 

710 FORMAT <6X, 15A4) 

NTOT =NF 

WRITE ( 1, 7£0) NTOT 
720 FORMAT (IX, 15) 

LEN ( 1 ) =X ( 1 ) 

DO 730 N=£, NTOT 
730 LEN (N) =X (N) —X (N— 1 ) 

WRITE (1,740) (LEN(N),EI (N) , WT (N) , N=1 , NTOT) 

740 FORMAT ( 1 X, F10. 6, F20. 1,7X,F15.6) 

WRITE (1, 750) EIO 
750 FORMAT ( 1 X, F20. 1) 

WRITE (1, 760) RRHUB 
760 FORMAT ( IX, F20. 3) 

C 

C CHECK THE TOTAL WEIGHT OF THE DISCRETE MODEL AND THE LOCATION 

C OF THE CG OF THE MODEL WITH RESPECT TO THE BASE FLANGE OF THE 

C SPOOL PIECE ( X =0 ) . 
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c 


c 


c 


c 


aoo 


aio 

aao 

830 

840 


850 


860 


870 

880 

890 


WT0T=0. 0 
WX=0. 0 
WSP=0. O 
WBD=0. O 
WTT=0. O 
WXSP=0. O 
WXBD=0. O 
WXTT=0.0 
DO 800 N=1 , NTOT 
WTOT=WTOT+WT(N> 

IF(X (N) . LE. 17. 75) THEN 

SPOOL PIECE SECTION 
WXSP=WXSP+WT <N) *X (N) 

WSP=WSP+WT(N) 

ELSE IF (X (N) . LE. 459. O) THEN 
BLADE SECTION 
WXBD=WXBD+WT (N> *X (N) 

WBD=WBD+WT(N> 

ELSE 

TRANSITION AND TIP SECTIONS 
WXTT=WXTT+WT (N) *X (N) 

UJTT=WTT+WT<N) 

END IF 

WX=WX+WT<N)*X(N) 

XB=WX/WTOT 
XSP=WXSP/WSP 
XBD=WXBD/WBD-17. 75 
XTT=WXTT /WTT— 459. 0 
WRITE (*, 810) 

FORMAT (//, SX, 1SHT0TAL WEIGHT, 1SX, 1 5HCG FROM HUB END) 

WRITE (*, 820) 

FORMAT (5X, 5H<LBS> , SOX, 6H<INCH) ) 

WRITE <*, 830) WTOT, XB 
FORMAT <3X, F9. 3, 16X, F8. 3) 

WRITE (*, 840) 

FORMAT (//, 2X, 1SHSP00L WEIGHT, 1SX, 15HCG FROM HUB END) 

WRITE <*, 820) 

WRITE (*, 830) WSP, XSP 
WRITE (*, 850) 

FORMAT (//, 2X, 12HBLADE WEIGHT, 12X, 17HCG FROM SPOOL END) 

WRITE <*, 820) 

WRITE <*, 830) WBD, XBD 
WRITE (*, 860) 

FORMAT <//, 2X, 19HTRANS. & TIP WEIGHT, 5X, 17HCG FROM BLADE END) 
WRITE <*, 820) 

WRITE (*, 830) WTT, XTT 
WRITE (*,870) 

FORMAT (//, 2X, 36HT0TAL MOMENT ABOUT CENTERLINE OF HUB) 

WRITE <*, 880) 

FORMAT < 15X, lOH < LBS- X NCH) > 

WRITE (*,890) WTOT* ( XB+22. 0) 

FORMAT (15X, F10. 1) 

STOP 

END 
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$nof loatcal Is 


^storage : £ 


C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


A. 2 CWTFLAP - counterweight assembly model generation code 

THIS PROGRAM GENERATES THE INPUT DATA FOR THE SPOOL 
PIECE (CF 764£54), THE TAPERED STEEL SPAR EXTENSION 
(CF 760549) AND THE COUNTERWEIGHT (CF 764554) ASSEMBLY FOR 
FLAPWISE VIBRATION (VIBRATION PERPENDICULAR TO THE PLANE OF 
ROTATION). THE PROGRAM CONSIDERS THE TAPERED SPAR EXTENSION 
TO BE BROKEN UP INTO 6 SECTIONS AS FOLLOWS : 


SECTION 1 - THE FLANGE AT THE BASE OF THE SPOOL PIECE 

OF LENGTH £. £5" WEIGHING APPROX. 108.74 LBS. 

SECTION £ - THE WEBBED MAIN BODY OF THE SPOOL PIECE 

OF LENGTH 13.75" WEIGHING APPROX. £07.96 LBS. 

SECTION 3 - THE FLANGE AT THE SPAR END OF THE SPOOL PIECE 
OF LENGTH 1.75" WEIGHING APPROX. 8£. 78 LBS. 

SECTION 4 - THE FLANGE AT THE BASE OF THE SPAR EXTENSION 
OF LENGTH 1.80" WEIGHING APPROX. 57.86 LBS. 

SECTION 5 - THE TRANSITION SECTION BETWEEN THE BASE 

FLANGE AND THE TAPERED PORTION OF THE SPAR. 
THIS SECTION IS 5. £0" IN LENGTH AND WEIGHS 
56. 63 LBS. 


SECTION 6 - THE TAPERED PORTION OF THE SPAR EXTENSION 
OF LENGTH 103.80" WEIGHING APPROX. 

795.96 LBS. 

SECTION 7 - THE TRANSITION SECTION BETWEEN THE TAPERED 

PORTION OF THE SPAR EXTENSION AND THE FLANGE 
AT THE TIP OF THE SPAR. THIS SECTION IS 
3.40" IN LENGTH AND WEIGHS 46.13 LBS. 


SECTION 8 - THE FLANGE AT THE TIP OF THE SPAR EXTENSION. 

THIS SECTION IS 1.80" IN LENGTH AND WEIGHS 
1£1. 56 LBS. 


CW THE LUMPED WEIGHT OF THE COUNTERWEIGHT AT THE FREE END OF 

THE STEEL SPAR EXTENSION (LBS) 

El (N) THE PRODUCT OF MODULUS OF ELASTICITY AND AREA MOMENT 
OF INERTIA AT MASS N (LBS*INCH*INCH) 

EIO THE PRODUCT OF MODULUS OF ELASTICITY AND AREA MOMENT 

OF INERTIA AT ZERO OR THE FIXED END (LBS*INCH*INCH> 
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c 

p 

L (N) 

THE 

LENGTH OF SECTION TO LEFT 

OF MASS 

N (INCHES) 

w 

c 

p 

NB 

THE 

MASS NUMBER AT THE BEGINNING OF A 

SECTION 

L# 

c 

p 

NF 

THE 

MASS NUMBER AT THE END OF 

A SECTION 

L 

c 

p 

ND 

THE 

NUMBER OF DIVISIONS IN A 

PARTICULAR SECTION 

c 

p 

NS 

THE 

NUMBER OF THE SECTION UNDER CONSIDERATION 

L, 

c 

P 

NTOT 

THE 

TOTAL NUMBER OF MASSES 



L 

c 

P 

RRHUB 

THE 

RADIUS OF THE RIGID HUB (INCHES) 


L/ 

c 

P 

TITLE 

THE 

PROBLEM TITLE STATEMENT 



L/ 

c 

P 

WT (N) 

THE 

WEIGHT OF MASS N (LBS) 



L 

c 

X !N) 

THE 

AXIAL LOCATION OF MASS N 

(INCH) 



c 

c 

c 

IMPLICIT REAL*8(A-H, 0-Z) 

REAL*8 L, M, Ml 

DIMENSION El (700) , L <700) , TITLE < 15) , WT <700) , X (700) 
ASSIGN LOGICAL UNIT 1 TO ’VIBIN.DAT’ 

0PEN<1, FILE®’ VIBIN. DAT’ , STATUS®’ NEW’ ) 

READ DATA FROM TERMINAL 


5 


10 


£0 


30 


WRITE <*, 5) 

FORMAT <//, £X, ’ THIS PROGRAM GENERATES THE DATA FOR THE ’, 

1 ’ FLAPWISE ’ , /, £X, ’ VIBRATION OF THE COUNTERWEIGHT ASSEMBLY ’, 

£ ///, £X, ’ ENTER THE RADIUS OF THE RIGID HUB (INCHES) : ’ , *) 

READ ( *, *) RRHUB 
WRITE <*, 10) 

FORMAT!//’ SECTION 1 - SPOOL PIECE FLANGE L®2. £5" W=108. 74 LB’, 
1 //’ SECTION £ - SPOOL PIECE BODY L=13. 75" W=£07. 96 LB’, 

£ //’ SECTION 3 - SPOOL PIECE FLANGE L=1.75" W=8£. 78 LB’, 

3 //’ SECTION 4 - SPAR BASE FLANGE L=1.80" W=57. 86 LB’, 


4 

//’ 

SECTION 

5 

- SPAR 

TRANSITION 

PORTION 

L=5. £0" W=5£. 63 

L-B’ , 

5 

//’ 

SECTION 

6 

- SPAR 

TAPERED PORTION L= 

103.80" W=795. 96 

LB’ , 

6 

//’ 

SECTION 

7 

- SPAR 

TRANSITION 

PORTION 

L=3. 40" W=46. 13 

LB’ , 

7 

//’ 

SECTION 

8 

- SPAR 

TIP FLANGE 

L=l. 80" 

W=l£ 1 . 56 LB’ ) 



WRITE!*, £0) 

FORMAT!/, 4X, ’ THE MAXIMUM TOTAL NUMBER OF DIVISIONS YOU ’ , 

1 / , 4X, ’ ARE ALLOWED IN ALL 8 SECTIONS IS 699. ’ ) 

£=30000000. O 
PI=3. 141 59S654 
NS® 1 
NF=0 

GO TO <40, 60, 80, 100, £00, 300, 400, 500, 600) NS 
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c 

C SECTION 1 PARAMETERS 

C 

40 WRITE (*,50) 

50 FORMAT (//, » ENTER THE NUMBER OF DIVISIONS FOR SECTION 1 : ’,$) 

READ ( *, 120)ND 
WC=48. 3268 
XL=0. 

XR=2. 25 
D0=21. 63 
DI=15. 825 
SO TO 520 
C 

C SECTION 2 PARAMETERS 

C 

60 WRITE!*, 70) 

70 FORMAT < //, • ENTER THE NUMBER OF DIVISIONS FOR SECTION 2 : » , *) 

READ ( *, 120) ND 
WC=15. 1242 
XR=16. 0 
C 

C CALCULATIONS FOR SECTION 2 

C 

DN=ND 

DX= ( XR— XL) /DN 
XR=XL+DX 
NB=NF+1 
NF=NF+ND 
SEI=E*2106. 268 
DO 75 N=NB, NF 
WT <N) =WC*DX 
X (N)=. 50* (XR+XL) 

El (N) =SEI 
XL=XR 

75 XR=XL+DX 

NS=NS+1 
GO TO 30 
C 

C SECTION 3 PARAMETERS 

C 

80 WRITE (*, SO) 

90 FORMAT!//, ’ ENTER THE NUMBER OF DIVISIONS FOR SECTION 3 : ’,*) 

READ! *, 120) ND 
WC=47. 3021 
XR=17. 75 
D0=21. 63 
DI=15. 97 
GO TO 520 
C 

C SECTION 4 PARAMETERS 

C 

100 WRITE!*, 110) 

110 FORMAT!//, * ENTER THE NUMBER OF DIVISIONS FOR SECTION 4 : ’,$) 

READ! *, 120) ND 
120 FORMAT! 15) 

WC=32. 14160448 
XR=19. 55 
DQ=20. 0 
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DI=15. 981 
GO TO 520 
C 

C SECTION 5 PARAMETERS 

C 

£00 WRITE (*,£10) 

£10 FORMAT (/, ’ ENTER THE NUMBER OF DIVISIONS FOR SECTION 5 : » , $) 

READ ( *, 1£0)ND 
M=-. 13980769£ 

B=£0. 46024038 
DI=16. O 
C0=. £83*PI/4. 

C1=C0*M**2 
C£=C0*2. *M*B 
C3=C0* ( B**£-D 1**2 ) 

XR=£4. 75 
GO TO 4£0 
C 

C SECTION 6 PARAMETERS 

C 

300 WRITE (*,310) 

310 FORMAT </, ’ ENTER THE NUMBER OF DIVISIONS FOR SECTION 6 : » , *) 

READ ( *, 1£0)ND 
M=. 014450867 
B=16. 64£34105 
T~. 50 

Ml=. £83*PI*T*M 
Bl = . £83*PI*T* (B— T) 

XR=1£8. 55 
C 

C CALCULATIONS FOR SECTION 6 

C 

DN=ND 

DX= ( XR— XL) /DN 

XR=XL+DX 

NB=NF+1 

NF=NF+ND 

DO 3£0 N=NB, NF 

WT (N) = (XR**£-XL**£> *M1 /£. O+Bl* < XR-XL) 

X (N)=( ( XR**3— XL**3) *M 1/3. 0+ < XR**£-XL**£) *B1 /£. O) /WT(N) 

DO=M*X (N)+B 

El (N) = <D0**3-3. 0*T*D0**£+4. 0*D0*T**2-£. 0*T**3> *PI*T*E/8. 0 
XL=XR 

3£0 XR=XL+DX 

NS=NS+1 
GO TO 30 
C 

C SECTION 7 PARAMETERS 

C 

400 WRITE (*,410) 

410 FORMAT ( /, ’ ENTER THE NUMBER OF DIVISIONS FOR SECTION 7 : ’ , *) 

READ ( *, 1£0>ND 
M=. 388823530 
B=-31. 48326478 
DI = 17.5 
C1=C0*M**2 
C£=CO*£. *M*B 
C3=C0* (B**2— DI**£) 
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XR=131. 95 
C 

c CALCULATIONS FOR SECTIONS 5 OR 7 

C 

420 DN=ND 

DX= ( XR— XL) /DN 

XR=XL+DX 

NB=NF+1 

NF=NF+ND 

DO 430 N=NB, NF 

WT <N) = ( X R**3— XL**3) *Cl/3. 0+ ( X R**£-XL**£) *C2/2. 0+C3* ( XR-XL) 

X <N) = ( (XR**4-XL**4) *C 1/4. 0+ ( XR**3-XL**3) *C2/3. 0 
1 +(XR**2-XL**2)*C3/£. 0)/WT(N) 

DO=M*X <N)+B 

El (N) - (E*PI* (D0**4-DI**4> ) /64. 0 
XL=X R 

430 XR=XL+DX 

NS=NS+1 
GO TO 30 
C 

C SECTION 8 PARAMETERS 

C 

500 WRITE (*,510) 

510 FORMAT (/, ’ ENTER THE NUMBER OF DIVISIONS FOR SECTION 8 : ’ $) 

READ ( *, 120) ND 
WC=67. 53381198 
XR=133. 75 
D0=24. 7 
DI = 1 7. 5 
C 

C CALCULATIONS FOR SECTIONS 1, 3, 4 OR 8 

C 

520 DN=ND 

DX= < XR— XL) /DN 
XR=XL+DX 
NB=NF+1 
NF=NF+ND 

SEI= (E*PI* (D0**4— DI**4) ) /64. 0 
DQ 530 N=NB, NF 
WT (N) =WC* (XR— XL) 

X (N)=. 50* (XR+XL) 

El (N)=SEI 
XL=XR 

530 XR=XL+DX 

NS=NS+1 
GO TO 30 

ACCOUNT FOR THE LUMPED COUNTERWEIGHT AT THE FREE END OF THE 
SPAR EXTENSION AND CALCULATE El AT THE BASE OR FIXED END OF 
THE SPOOL PIECE (X=0. O). 

600 DQ=£1 . 63 

D I = 1 5. 8£5 

EIO= (E*PI* (D0**4-DI**4) ) /64. O 
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WRITE(*, 610) 

610 FORMAT ( /, * ENTER THE LUMPED WEIGHT OF THE COUNTERWEIGHT AT THE 

1 /, * FREE END OF THE TAPERED SPAR EXTENSION <LBS>: * , %) 

READ (*, *) CW 
WRITE <*, 630) 

630 FORMAT ( /, ’ ENTER A TITLE STATEMENT (60 CHAR. MAX.) : ’ , $) 

READ (*, 640) TITLE 
640 FORMAT (15A4) 

C IF NO LUMPED COUNTERWEIGHT EXISTS, OUTPUT THE DATA FOR THE TAPEREI 

C STEEL SPAR MODEL. IF A LUMPED COUNTERWEIGHT EXISTS, ADD ONE MORE 

C LUMPED WEIGHT TO THE MODEL AT THE CG OF THE COUNTERWEIGHT, THEN 

C ADD AN ADDITIONAL BEAM SEGMENT TO THE MODEL. 

C 

IF (CW. EQ. 0. O) GO TO 700 
C 

C ACCOUNT FOR THE COUNTERWEIGHT AT THE END OF THE SPAR EXTENSION 

C 

NF=NF+1 
WT(NF)=CW 
A=19. 34 
B= 1 4. 50 

El (NF) = (E*PI*A*B**3> /4. O 
X (NF) =133. 75+9. 8029 
C 

C OUTPUT RESULTS TO VIBIN DATA FILE 

C 

700 WRITE (1, 710) TITLE 

710 FORMAT (6X, 15A4) 

NTOT=NF 

WRITE (1, 720) NTOT 
720 FORMAT (IX, 15) 

L ( 1 > =X ( 1 ) 

DO 730 N=2, NTOT 
730 L(N)=X(N>— X(N-l) 

WRITE ( 1,740) (L (N) , El (N) , WT (N) , N=1 , NTOT) 

740 FORMAT (IX, F10. 6, F20. 1, 7X, F15. 6) 

WRITE (1, 750) EIO 
750 FORMAT ( 1 X, F20. 1) 

WRITE (1, 760) RRHUB 
760 FORMAT (IX, F20. 3) 

C 

C CHECK THE TOTAL WEIGHT OF THE DISCRETE MODEL AND THE LOCATION 

C OF THE CG OF THE -MODEL WITH RESPECT TO THE BASE FLANGE OF THE 

C SPOOL PIECE (X=0). 

C 

WT0T=0. 0 
WX=0. 0 

DO 800 N= 1 , NTOT 
WTOT=WTOT+WT (N) 

800 WX=WX+WT(N)*X (N) 

XB=WX/WTOT 
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810 

820 

830 

840 

850 

860 


WRI TE (*, 810) 

FORMAT (/, 2X, 12HT0TAL WEIGHT, 5X, 2HCG) 

WRITE (■*, 820) 

FORMAT (5X,5H (LBS), 7X,6H( INCH) ) 

WRITE (■*, 830) WTOT, XB 
FORMAT (3X, F3. 3, 3X, F8. 3) 

WRITE <■*, 840) 

FORMAT (//, 2X, 36HT0TAL MOMENT ABOUT CENTERLINE OF HUB) 
WRITE (•*, 850) 

FORMAT ( 15X, 10H (LBS-INCH) ) 

WRITE (*, 860) WTOT* (XB+22. 0) 

FORMAT (15X, F10. 1) 

STOP 

END 
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* 


* 


A. 3 CWTCHRD counterweight assembly model generation code 

THIS PROGRAM GENERATES THE INPUT DATA FOR THE SPOOL 
PIECE (CF 764254), THE TAPERED STEEL SPAR EXTENSION 
<CF 760549) AND THE COUNTERWEIGHT (CF 764554) ASSEMBLY FOR 
CHORDWISE VIBRATION (VIBRATION IN THE PLANE OF 


*NOFLOATCALLS 
^•STORAGE : 2 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


ROTATION). THE 
TO BE BROKEN UP 


PROGRAM CONSIDERS THE TAPERED SPAR EXTENSION 
INTO Q SECTIONS AS FOLLOWS : 


SECTION 1 - THE FLANGE AT THE BASE OF THE SPOOL PIECE 

OF LENGTH 2.25" WEIGHING APPROX. 108.74 LBS. 

SECTION 2 - THE WEBBED MAIN BODY OF THE SPOOL PIECE 

OF LENGTH 13.75" WEIGHING APPROX. 207.96 LBS. 

SECTION 3 - THE FLANGE AT THE SPAR END OF THE SPOOL PIECE 
OF LENGTH 1.75" WEIGHING APPROX. 82.78 LBS. 

SECTION 4 - THE FLANGE AT THE BASE OF THE SPAR EXTENSION 
OF LENGTH 1.80" WEIGHING APPROX. 57.86 LBS. 


SECTION 


THE TRANSITION SECTION BETWEEN THE BASE 
FLANGE AND THE TAPERED PORTION OF THE SPAR. 
THIS SECTION IS 5.20" IN LENGTH AND WEIGHS 
56.63 LBS. 


SECTION 6 - THE TAPERED PORTION OF THE SPAR EXTENSION 
OF LENGTH 103.80" WEIGHING APPROX. 

795. 96 LBS. 


C 


C 

C 

C 

C 

C 

C 

C 


C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 


SECTION 7 - THE TRANSITION SECTION BETWEEN THE TAPERED 

PORTION OF THE SPAR EXTENSION AND THE FLANGE 
AT THE TIP OF THE SPAR. THIS SECTION IS 
3.40" IN LENGTH AND WEIGHS 46.13 LBS. 

SECTION 8 - THE FLANGE AT THE TIP OF THE SPAR EXTENSION. 

THIS SECTION IS 1.80" IN LENGTH AND WEIGHS 
121.56 LBS. 


CW THE LUMPED WEIGHT OF THE COUNTERWEIGHT AT THE FREE END OF 

THE STEEL SPAR EXTENSION (LBS) 

EI(N) THE PRODUCT OF MODULUS OF ELASTICITY AND AREA MOMENT 
OF INERTIA AT MASS N (LBS*INCH*INCH> 

EIO THE PRODUCT OF MODULUS OF ELASTICITY AND AREA MOMENT 

OF INERTIA AT ZERO OR THE FIXED END (LBS*INCH*INCH> 

L (N) THE LENGTH OF SECTION TO LEFT OF MASS N (INCHES) 
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c 

c 

NB 

THE 

c 

c 

NF 

THE 

c 

c 

ND 

THE 

c 

c 

NS 

THE 

c 

c 

NTOT 

THE 

c 

c 

RRHUB 

THE 

c 

c 

TITLE 

THE 

c 

c 

WT (N) 

THE 

c 

X (N) 

THE 


c 

c 

c 


MASS NUMBER AT THE BEGINNING OF FI SECTION 
MASS NUMBER AT THE END OF A SECTION 
NUMBER OF DIVISIONS IN A PARTICULAR SECTION 
NUMBER OF THE SECTION UNDER CONSIDERATION 
TOTAL NUMBER OF MASSES 
RADIUS OF THE RIGID HUB (INCHES) 

PROBLEM TITLE STATEMENT 
HEIGHT OF MASS N (LBS) 

AXIAL LOCATION OF MASS N (INCH) 


IMPLICIT R£AL*8 ( A— H, 0— Z ) 

REAL*8 L, M, Ml 

DIMENSION El (700) , L(700) , TITLE (15) , WT(700> , X (700) 
ASSIGN LOGICAL UNIT 1 TO ’VIBIN.DAT’ 


10 


£0 


1 

£ 


OPEN ( 1 , FILE = ’ VIBIN. DAT’ , STATUS=’ NEW’ ) 

READ DATA FROM TERMINAL 
WRITE!*, 5) 

FORMAT (//, EX, ’ THIS PROGRAM GENERATES THE DATA FOR THE 
’CHORDWISE ’ , /, £X, ’ VIBRATION OF THE COUNTERWEIGHT ASSEMBLY ’, 
// /, £X, ’ ENTER THE RADIUS OF THE RIGID HUB (INCHES) : ’ , *> 

READ (*. *) RRHUB 
WRITE (*, 10) 

FORMAT (//2X, ’ SECTION 1-SPOOL PIECE FLANGE L=2. £5" W=108. 74 LB’, 


//2X, 

//£X, 

//£X, 

//£X, 

//£X, 

//£X, 

//£X 


SECTION E-SPOOL PIECE BODY L=13. 75" W=207. 96 LB’, 

SECTION 3-SPOOL PIECE FLANGE L=1.75" W=8£. 78 LB’, 

SECTION 4-SPAR BASE FLANGE L=1.80" W=57. 88 LB’, 

SECTION 5-SPAR TRANSITION PORTION L=5. £0" W=5£. 63 LB’, 


30 


SECTION 7-SPAR TRANSITION PORTION L=3. 40" W=46. 13 LB’, 

SECTION 8-SPAR TIP FLANGE L=1.80" W=l£1.56 LB’) 

WRITE!*, £0) 

FORMAT!/, 4X, ’ THE MAXIMUM TOTAL NUMBER OF DIVISIONS YOU ’ , 

1 /, 4X, ’ ARE ALLOWED IN ALL 8 SECTIONS IS 699. ’ ) 

E=30000000. O 
PI=3. 141598654 
NS=1 
NF=0 

GO TO (40, 60, 80, 100, £00, 300, 400, 500, 600) NS 
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c 

C SECTION 1 PARAMETERS 

C 

AO WRITE <*, 50) 

50 FORMAT <//, ’ ENTER THE NUMBER OF DIVISIONS FOR SECTION 1 : ’ , *) 

READ <*, 120) ND 
WC=A8. 3266 
XL=0. 

XR=2. 25 
D0=21. 63 
DI— 15. 825 
SO TO 520 
C 

C SECTION 2 PARAMETERS 

C 

60 WRITE (*,70) 

70 FORMAT <//, ’ ENTER THE NUMBER OF DIVISIONS FOR SECTION 2 : ’ , $> 

READ (*, 120) ND 
WC=15. 1242 
XR=16. 0 
C 

C CALCULATIONS FOR SECTION 2 

C 

DN=ND 

DX= ( XR-XL) /DN 
XR=XL+DX 
NB=NF+1 
NF=NF+ND 
SEI=E*2106. 268 
DO 75 N=NB, NF 
WT(N)=WC*DX 
X (N)=. 50* (XR+XL) 

El <N)=SEI 
XL=XR 

75 XR=XL+DX 

NS=NS+i 
GO TO 30 
C 

C SECTION 3 PARAMETERS 

C 

80 WRITE (*, 90) 

90 FORMAT ( //, ’ ENTER THE NUMBER OF DIVISIONS FOR SECTION 3 : ’,*) 

READ (*, 120) ND 
WC=47. 3021 
XR=17. 75 
D0=21. 63 
D I = 1 5. 97 
GO TO 520 
C 

C SECTION A PARAMETERS 

C 

100 WRITE (*,110) 

110 FORMAT (//, ’ ENTER THE NUMBER OF DIVISIONS FOR SECTION A : ’ , *> 

READ (*, 120) ND 
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120 


C 

c 

c 

£00 

£10 


c 

CJ 

c 

300 

310 


C 

c 

c 


320 


C 

r* 

c 

A00 

A10 


FORMAT (15) 

WC=3£. 1 A160AAQ 
XR=19. 55 
D0=20. 0 
DI=15. 981 
GO TO 520 

SECTION 5 PARAMETERS 
WRITE (*, £10) 

FORMAT (/, » ENTER THE NUMBER OF DIVISIONS FOR SECTION 5 ; ’,*) 

READ (*, 120) ND 

M=— . 139807692 

B=£0. A602A038 

DI=16. 0 

C0=. £S3*PI/A. 

C1=C0*M**2 
C2=C0*£. *M*B 
03=00* (B**£-DI**£) 

XR=£A. 75 
GO TO A20 

SECTION 6 PARAMETERS 
WRITE (*, 310) 

FORMAT </, ’ ENTER THE NUMBER OF DIVISIONS FOR SECTION 6 : ’ , *> 

READ (*, 120) ND 

M=. Ol AA50867 

B=16. 6A23A105 

T=. 50 

Ml=. £83*PI*T*M 
Bl=. 283*PI*T*(B-T) 

XR— 128. 55 


CALCULATIONS FOR SECTION 8 
DN=ND 

DX= < XR— XL) /DN 

XR=XL+DX 

NB=NF+1 

NF=NF+ND 

DO 320 N=NB, NF 

WT <N) = (XR**2-XL**2) *M1 /2. O+Bl* ( XR-XL) 

X (N) =< (XR**3-XL**3) *Ml/3. 0+ ( XR**£-XL**2) *Bl/2. 0) /WT (N) 
DO=M*X(N)+B 

El (N) = (D0**3— 3. 0*T*D0**2+A. 0*D0*T**2-£. 0*T**3) *PI*T*E/8. O 

XL=XR 

XR=XL+DX 

NS=NS+1 

GO TO 30 

SECTION 7 PARAMETERS 
WRITE (*, A10) 

FORMAT (/, ’ ENTER THE NUMBER OF DIVISIONS FOR SECTION 7 : ’,*) 

READ (*, 120) ND 
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M=. 388823530 
B=-31. 48326478 
1)1 = 17. 5 
C1=C0*M**£ 

C£=C0*£. *M*B 
C3=C0* (B**2-DI**2> 

XR=131. 95 
C 

C CALCULATIONS FOR SECTIONS 5 OR 7 

C 

420 DN=ND 

DX= ( XR— XL) /DN 

XR=XL+DX 

NB=NF+1 

NF=NF+ND 

DO 430 N=NB, NF 

WT (N) =(XR**3-XL**3)*Cl/3. 0+ ( XR**£-XL**£) *C£/£. 0+C3* (XR-XL) 

X (N) = < < XR**4— XL**4) *Cl/4. 0+ ( XR**3-XL**3) *C£/3. 0 
1 + (XR**£-Xi_**£) *C3/£. O) /WT (N) 

DO=M*X (N> +B 

El (N) = (E*PI* (D0**4-DI**4) ) /64. O 
XL=XR 

430 XR=XL+DX 

NS=NS+1 
SO TO 30 
C 

L SECTION 8 PARAMETERS 

C 

500 WRITE (*,510) 

510 FORMAT ( /, ’ ENTER THE NUMBER OF DIVISIONS FOR SECTION 8 : ’,*) 

READ (*, 120) ND 
WC=S7. 53381198 
X R= 1 33. 75 
D0=£4. 7 
D I = 1 7. 5 
C 

C CALCULATIONS FOR SECTIONS 1, 3, 4 OR 8 

C 

520 DN=ND 

DX= (XR— XL) /DN 
X R=XL+DX 
NB=NF+ 1 
NF=NF+ND 

SEI= (E*PI* (D0**4— DI**4) ) /64. 0 
DO 530 N=NB, NF 
WT(N)=WC* (XR-XL) 

X (N>=. 50* ( XR+XL) 

El (N)=SEI 
XL=XR 

530 XR=XL+DX 

NS=NS+1 
GO TO 30 
C 

C ACCOUNT FOR THE LUMPED COUNTERWEIGHT AT THE FREE END OF THE 

C SPAR EXTENSION AND CALCULATE El AT THE BASE OR FIXED END OF 

C THE SPOOL PIECE (X=0. O). 
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600 D0=£1.63 

DI= 1 5. 82 5 

EI0= <E*PI* (D0**4-DI**4) ) /64. 0 
WRITE (*, 610) 

610 FORMAT (/, * ENTER THE LUMPED WEIGHT OF THE COUNTERWEIGHT AT THE 

1 /, ’ FREE END OF THE TAPERED SPAR EXTENSION <LBS) : ’ ,*) 

READ <*, *)CW 
WRITE <*, 630) 

630 FORMAT (/, ’ ENTER A TITLE STATEMENT (60 CHAR. MAX.) : » , *) 

READ (*, 640) TITLE % 

640 FORMAT (15A4) 

C IF NO LUMPED COUNTERWEIGHT EXISTS, OUTPUT THE DATA FOR THE TAPEREE 

C STEEL SPAR MODEL. IF A LUMPED COUNTERWEIGHT EXISTS. ADD ONE MORE 

C LUMPED WEIGHT TO THE MODEL AT THE CG OF THE COUNTERWEIGHT, THEN * 

C ADD AN ADDITIONAL BEAM SEGMENT TO THE MODEL. 

C 

IF (CW. EQ. O. 0) GO TO 700 
C 

C ACCOUNT FOR THE COUNTERWEIGHT AT THE END OF THE SPAR EXTENSION 

C 

NF=NF+1 
WT <NF) =CW 
A= 1 9. 34 
B= 14. 50 

El <NF> = <E*PI*B*A**3> /4. 0 
X (NF) =133. 75+9. 8029 
C 

C OUTPUT RESULTS TO VIBIN DATA FILE 

C 

700 WRITE (1, 710) TITLE 

710 FORMAT (6X, 15A4) 

NTOT=NF 

WRITE (1, 720) NTOT 
720 FORMAT (IX, 15) 

L ( 1 ) =X ( 1 ) 

DO 730 N=£, NTOT 
730 L (N) =X (N) -X (N-l ) 

WRITE (1,740) (L (N) , El (N) , WT (N) , N=1 , NTOT) 

740 FORMAT ( IX, F10. 6, F20. 1,7X,F15. 6) 

WRITEd, 750)EI0 
750 FORMAT ( IX, F20. 1) 

WRITE (1,760) RRHUB 
760 FORMAT (IX, F20. 3) 

C 

C CHECK THE TOTAL WEIGHT OF THE DISCRETE MODEL AND THE LOCATION 

C OF THE CG OF THE MODEL WITH RESPECT TO THE BASE FLANGE OF THE 

C SPOOL PIECE (X=0). 

C 
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800 

810 

820 

830 

840 

850 

860 


WT0T=0. 0 
WX=0. 0 

DO 800 N=1 , NTOT 
WTOT=WTOT+WT<N) 

WX=WX+WT<N)*X(N) 

XB=WX/WTOT 
WRITE (*, 810) 

FORMAT (/, 2X, 12HT0TAL WEIGHT, 5X, 2HCG) 

WRITE <*, 820) 

FORMAT <5X, 5H(LBS) , 7X, 6H<INCH) ) 

WRITE (*, 830) WTOT, XB 
FORMAT (3X, F9. 3, 3X, F8. 3) 

WRITE (*, 840) 

FORMAT <//, EX, 36HT0TAL MOMENT ABOUT CENTERLINE OF HUB) 
WRITEt*, 850) 

FORMAT < 15X, 10H (LBS-INSH) ) 

WRITE (*, 860) WTOT* ( XB+22. O) 

FORMAT (15X, F10. 1) 

STOP 

END 
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$nof loatcal Is 
^Storage^ 

C 

C A. 4 BEAM - main analysis code 

C 

c 

C THIS IS THE MAIN PROGRAM USED FOR THE DETERMINATION OF THE NATURAL 

C FREQUENCIES AND MODE SHAPES FOR THE LATERAL VIBRATION OF A CANTILEVERED 

C BEAM. THE PROGRAM IS BASED ON A LUMPED MASS APPROXIMATION OF A 

C CONTINUOUS SYSTEM. THE LUMPED MASSES IN THE DISCRETE MODEL OF THE BEAM 

C ARE CONNECTED BY MASSLESS, BUT FLEXIBLE, BEAM SEGMENTS! THE FLEXURAL 

C RIGIDITY OF THE BEAM, El, CAN EITHER BE UNIFORM OR IT CAN VARY 

C ALONG THE LENGTH OF THE BEAM'. IF THE FLEXURAL RIGIDITY VARIES ALONG 

C THE LENGTH OF THE BEAM, A LINEAR VARIATION IN El IS ASSUMED IN 

C EACH BEAM SEGMENT BETWEEN TWO SUCCESSIVE MASSES. THE EFFECTS OF 

C ROTATION ON THE VIBRATORY CHARACTERISTICS OF THE BEAM CAN ALSO BE 

C DETERMINED FOR FLAPWISE VIBRATION USING THIS PROGRAM. 

C 

C 

C 

C 

C 

C El (N) THE PRODUCT OF MODULUS OF ELASTICITY AND AREA MOMENT 

C OF INERTIA AT MASS N <LBS*INCH*INCH> 

C 

C EIO THE PRODUCT OF MODULUS OF ELASTICITY AND AREA MOMENT 

C OF INERTIA AT ZERO OR THE FIXED END (LBS*INCH*INCH> 

C 

C ISPEED AN INDEX CORRESPONDING TO THE PARTICULAR ROTOR 

C SPEED UNDER INVESTIGATION 

C 

C L <N) THE LENGTH OF SECTION TO LEFT OF MASS N (INCHES) 

C 

C NPLACE THE NUMBER OF DIGITS REQUIRED TO THE RIGHT OF THE 

C DECIMAL POINT 

C 

C NRPM THE TOTAL NUMBER OF ROTATIONAL SPEEDS CONSIDERED 

C 

C N'fOT THE TOTAL NUMBER OF MASSES 

C 

C NW THE NUMBER OF NATURAL FREQUENCIES IN THE RANGE FROM 

C WD TO WL 

C 

C RPM(N) THE ROTOR ROTATIONAL SPEED ( REV/MIN) 

C 

C RRHUB THE RADIUS OF THE RIGID HUB (INCHES) 

C 

C TIME THE ELAPSED TIME DURING EXECUTION (SEC) 

C 

C TITLE THE PROBLEM TITLE STATEMENT 

C 

C W(I> THE NATURAL FREQUENCIES OF VIBRATION (RAD/SEC) 

C 

C WD THE INITIAL DELTA USED IN SEARCHING FOR NATURAL 

C FREQUENCIES (RAD/SEC) 
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c 

C WF THE FIRST VALUE TO BE CONSIDERED AS A POSSIBLE NATURAL 

C FREQUENCY (RAD/SEC) 

C 

C WL THE LAST VALUE TO BE CONSIDERED AS A POSSIBLE NATURAL 

C FREQUENCY (RAD/SEC) 

C 

C WT(N) THE WEIGHT OF MASS N (LBS) 

C 

C 

C 

* C 

IMPLICIT REAL*Q ( A— H, O— Z ) 

REAL*S L 

' CHARACTER*! KEY 

COMMON El (700) , EIO, FL(700) , L(700) , NPLACE, NRPM, NTOT, NW, RPM(IO) 

1 , RRHUB, TITLE ( 15) , W ( 100) , WD, WF, WL, WT (700) , KEY 

READ INITIAL DATA FROM TERMINAL AND INPUT FILE 

CALL INPUT 

CONSIDER EACH OF THE ROTATIONAL SPEED INDIVIDUALLY 
DO 10 ISPEED=1, NRPM 

CALCULATE NATURAL FREQUENCIES 
CALL NFREQ ( ISPEED) 

OUTPUT NATURAL FREQUENCY RESULTS 
CALL OUTPUT (ISPEED) 

CALCULATE AND OUTPUT MODE SHAPES IF DESIRED 
CALL MODES (ISPEED) 

CONTINUE 

TERMINATE 

STOP 
END 
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$nof loatcal is A. 5 DET - 1 subroutine to evaluate determinant 

*storage:£ 

REAL*Q FUNCTION DET(WW) 

C 

C THIS FUNCTION EVALUATES THE DETERMINANT WHICH REFLECTS THE 

C BOUNDARY CONDITIONS OF THE PROBLEM. WHEN THIS DETERMINANT 

t: IS ZERO, THE CORRESPONDING FREQUENCY IS A NATURAL FREQUENCY 

C OF THE SYSTEM. 

C 

C A ( I , J) THE MATRIX WHICH RELATES THE DEFLECTION, SLOPE, 

C MOMENT AND SHEAR AT THE N-l MASS TO THOSE QUANTITIES 

C AT MASS N 

C 

C NTOT THE TOTAL NUMBER OF MASSES 

C 

C U(I,J) THE CONCATENATION OF CAD MATRICIES RELATING DEFLECTION, 

C SLOPE, MOMENT AND SHEAR AT THE FIXED END TO THOSE 

C QUANTITIES AT THE FREE END 

C 

C UNEW ( I , J) SAME AS U(I,J) 

C 

C WW THE FREQUENCY AT WHICH THE DETERMINANT IS TO BE COMPUTED 

C 

IMPLICIT REAL*Q ( A-H, O-Z ) 

REAL*8 L 
CHARACTER*1 KEY 

COMMON El (700) , EIO, FL(700) , L (700) , NPLACE, NRPM, NTOT, NW, RPM(10> 

1 , RRHUB, TITLE (15) , W(100) , WD, WF, WL, WT(700) , KEY 

R£AL*8 A (4, 4) , U (4, 4) , UNEW (4, 4) 

C 

C INITIALIZE CUD (SET IT EQUAL TO THE IDENTITY MATRIX) 

C 

DO 10 1=1,4 
DO 10 J= 1,4 
U( I. J ) =o. 

IF(I.EQ.J) U(1,J>=1. 

10 CONTINUE 

C 

C COMPUTE CUD = THE CONCATENATION OF THE CAD MATRICES FOR ALL MASSES 

C 

DO £0 N=1 , NTOT 

C FILL CAD FOR THE MASS N 

CALL FILLA(A, N, WW) 

C CONCATENATE CAD WITH THE EXISTING CUD 

CALL MULT (A, 4, 4, 4, 4, U, 4, 4, 4, UNEW, 4, 4) 

C SET CUD = CUNEWD 

DO 15 1=1,4 
DO 15 J= 1 , 4 

15 U < I , J ) =UNEW ( I , J ) 

£0 CONTINUE 

C 

C COMPUTE DETERMINANT WHICH REFLECTS BOUNDARY CONDITIONS 

C 

DET=U (3, 3) *U (4, 4) — U (3, 4) *U (4,3) 

RETURN 

END 
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$nof loatcal Is A. 6 FILLA - subroutine to fill the A matrix 

$storage:£ 


C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 

c: 

c 

c 


c 

c 

c 

c 

c 

c 

c 


% 


SUBROUTINE FILLA (A, N, WW) 


THIS SUBROUTINE FILLS THE MATRIX A(I,J). THIS MATRIX IS 
THE MATRIX WHICH RELATES THE DEFLECTION, SLOPE, MOMENT AND 
SHEAR AT THE N-l MASS TO THOSE QUANTITIES AT MASS N. 


A < I , J) THE MATRIX WHICH RELATES THE DEFLECTION, SLOPE, MOMENT 
AND SHEAR AT THE N-l MASS TO THOSE QUANTITIES AT MASS N 

AM THE SLOPE AT THE FREE END OF A CANTILEVERED BEAM 

DUE TO A UNIT MOMENT APPLIED AT THE FREE END (RAD/LBS-IN) 

AS THE SLOPE AT THE FREE END OF A CANTILEVERED BEAM 

DUE TO A UNIT FORCE APPLIED AT THE FREE END (RAD/LBS) 

DM THE DEFLECTION AT THE FREE END OF A CANTILEVERED BEAM 

DUE TO A UNIT MOMENT APPLIED AT THE FREE END (1/LBS) 

DS THE DEFLECTION AT THE FREE END OF A CANTILEVERED BEAM 

DUE TO A UNIT FORCE APPLIED AT THE FREE END (INCH/LBS) 

EIAVE THE AVERAGE VALUE OF El AT MASS N, EI(N), AND El AT 
MASS N-l, El (N-l) (LBS*- INCH* INCH) 

EIN THE VALUE OF THE FLEXURAL RIGIDITY AT MASS N, El (N) 

(LBS* INCH* INCH) 

EINM1 THE VALUE OF THE FLEXURAL RIGIDITY AT MASS N-l, El (N-l) 
(LBS* INCH* INCH) 

EIO THE VALUE OF THE FLEXURAL RIGIDITY AT ZERO OR THE FIXED 

LN THE DISTANCE BETWEEN MASS N-l AND MASS N (INCHES) 

MASSN THE MASS OF THE N TH LUMPED MASS (LBS-SEC*SEC/ IN) 

RPM THE ROTOR ROTATIONAL SPEED ( REV/MIN) 

RRHUB THE RADIUS OF THE RIGID HUB (INCHES) 

WW THE CURRENT FREQUENCY AT WHICH THE DETERMINANT IS TO BE 

COMPUTED (RAD/SEC) 


IMPLICIT REAL*8 ( A— H, O-Z ) 

REAL*0 A (A, A) , LN, MASSN, L 
CHARACTER* 1 KEY 

COMMON El (700) , EIO, FL (700) , L (700) , NPLACE, NRPM, NTOT, NW, RPM ( 10) 
, RRHUB, TITLE ( 15) , W ( 100) , WD, WF, WL, WT (700) , KEY 
IF (N. GT. 1) GO TO 10 
EINM1=EI0 
A ( 1 , 1 ) — 1 . O 
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ft (£, 1)=0.0 
0(3, 1 > =0. 0 
GO TO £0 

10 EINM1=EI (N— 1) 

£0 EIN=EI (N> 

MASSN=WT (N) /3Q6. 4 
LN=L(N> 

El AVE= (EIN+EINMi ) /£. 0 
C 

C IF THE VALUE OF El (N) /El (N-l ) IS BETWEEN . 9999 AND 1.00009 THEN 

THE BEAM IS CONSIDERED TO BE UNIFORM BETWEEN THE' TWO MOSSES WITH > 
ON OVERAGE VALUE OF El. IF NOT, O LINEAR VARIATION IN El BETWEEN 
THE TWO MOSSES IS ASSUMED. 

R=EIN/EINM1 

IF (R. GT. . 9999. AND. R. LT. 1. 00009) GO TO 30 

ASSUMING A LINEAR VARIATION IN El BETWEEN El (N) AND EI(N-l) 

DEI=EIN— EINM1 
REI=EIN/EINM1 

DS= < (LN**3) /DEI )*( (1. 0+£. 0*EINM1/DEI+ (EINM1 /DEI ) *•*£) * 

1 DLOG < RE I ) -1 . 5-EINM1 /DEI ) 

DM=( (LN**£) /DEI ) * ( (E IN/DE I ) *DLOG ( RE I ) -1 . O) 

AM= (LN/DEI ) *DLOG ( REI ) 

GO TO 40 

ASSUMING A UNIFORM VALUE OF El BETWEEN MASSES 

30 DS= (LN**3) / <3. 0*EIAVE> 

DM= (LN**£) / <£. 0*EI AVE) 

OM=LN/EIAVE 

C 

40 AS=DM 

DEN=1 . + (AS-DM) *FL (N) + (AM*DS-AS*DM) *FL <N) **£ 

0(4, 1 ) =WW**£*MASSN 
A C 1 , £ ) = ( LN- ( DS-LN* AS ) *FL < N ) ) /DEN 
A (£,£) = ( 1. - (DM-LN*OM) *FL (N) ) /DEN 
0(3, £>=A(1, £>*FL(N> 

0(4, £)=WW**£*M0SSN*A(1, £) 

0(1,3)= (DM- (AM*DS— AS*DM> *FL (N) ) /DEN 

A (£, 3) =AM/DEN 

0(3, 3)= ( 1. +AS*FL (N) ) /DEN 

A (4, 3) =WW**£*MASSN*A (1,3) 

0(1, 4) = (LN*DM— DS— LN* (AM*DS-AS*DM) *FL (N) ) /DEN 
A (£, 4) = <LN*OM-AS— <OM*DS-AS*DM> *FL (N) ) /DEN 
0(3, 4)=A(1, £) 

A (4, 4) =WW**£*MASSN*A ( 1 , 4) +1. 

RETURN 

END 
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$nof loatcai Is A. 7 FORCE -*• subroutine to calculate centrifugal forces 

^storage : 2 

SUBROUTINE FORCE ( ISPEED) 

C 

C 

C THIS SUBROUTINE CALCULATES THE CENTRIFUGAL FORCE ACTING TO THE 

C LEFT OF EACH MASS. 

C 
C 

L FL (N) 

C 
C 

C ISPEED 

C 
C 

C L <N) 

C 

C NTOT 

C 

C OMEGA 

C 

C RPM< ISPEED) 

C 

C RRHUB 

C 

C WT(N) 

C 

C XI 

c 
c 
c 
c 
c 

C COMPUTE THE LENGTH TO THE LAST MASS 

C 

IMPLICIT REAL*8 ( A— H, O— Z > 

REAL*8 L 
CHARACTERS KEY 

COMMON El (700) , EIO, FL(700) , L(700) , NPLACE, NRPM, NTOT, NW, RPM(IO) 

1 , RRHUB, TITLE ( 15), W( 100), WD, WF, WL, WT (700), KEY 

C 

XN=RRHUB 
DO 10 1=1, NTOT 
10 XN=XN+L(I) 

C 

C COMPUTE THE INERTIAL FORCES BEGINNING AT THE LAST MASS AND 

C PROCEEDING TO MASS 1 

C 
C 

OMEGA=RPM( ISPEED) *3. 14159265/30. 

X I =XN 

FL (NTOT) =WT (NTOT) /386. 4*XN*0MEGA**£ 

DO 20 1=2, NTOT 
J=NT0T+1— I 
X I=X I— L ( J+l ) 

20 FL ( J ) =FL ( J+l ) +WT ( J ) /386. 4*X I *OMEGA**£ 

RETURN 
END 


THE CENTRIFUGAL FORCE ACTING TO THE LEFT OF 
THE N TH MASS 

AN INDEX CORRESPONDING TO THE PARTICULAR ROTOR 
SPEED UNDER INVESTIGATION 

THE DISTANCE BETWEEN MASS N-l AND MASS N (INCHES) 

THE TOTAL NUMBER OF MASSES 

THE ROTOR ROTATIONAL SPEED (RAD/SEC) 

THE CURRENT ROTOR ROTATIONAL SPEED ( REV /MIN) 

THE RADIUS OF THE RIGID HUB (INCHES) 

THE WEIGHT OF MASS N (LBS) 

THE DISTANCE FROM THE CENTERLINE OF THE ROTOR 
TO THE N TH MASS (INCHES) 
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$nof loatcal Is 
istorage : £ 


A. 8 INPUT - subroutine to perform all input functions 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

CJ 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE INPUT 

THIS SUBROUTINE PERMITS THE USER TO ENTER INITIftL VALUES FROM 
THE TERMINAL AND ALSO READS BEAM DATA FROM THE FILE 
’ VIBIN. DAT’ . 


El (N) 

EIO 

KEY 

L(N) 

NPLACE 

NRPM 

NTOT 

RPM(N) 

RRHUB 

TITLE 

WD 

WF 

WL 

WT(N) 


THE PRODUCT OF MODULUS OF ELASTICITY AND AREA MOMENT 
OF INERTIA AT MASS N (LBS*INCH*INCH) 

THE PRODUCT OF MODULUS OF ELASTICITY AND AREA MOMENT 
OF INERTIA AT ZERO OR THE FIXED END <LBS*INCH*INCH> 

A FLAG WHICH IF “Y" WILL CAUSE THE MODE SHAPES TO BE 
CALCULATED 

THE LENGTH OF SECTION TO LEFT OF MASS N (INCHES) 

THE NUMBER OF DIGITS REQUIRED TO THE RIGHT OF THE 
DECIMAL POINT 

THE TOTAL NUMBER OF ROTATIONAL SPEEDS CONSIDERED 

THE TOTAL NUMBER OF MASSES 

THE ROTOR ROTATIONAL SPEED (REV/MIN) 

THE RADIUS OF THE RIGID HUB (INCHES) 

THE PROBLEM TITLE STATEMENT 

THE INITIAL DELTA USED IN SEARCHING FOR NATURAL 
FREQUENCIES (RAD/SEC) 

THE FIRST VALUE TO BE CONSIDERED AS A POSSIBLE NATURAL 
FREQUENCY (RAD/SEC) 

THE LAST VALUE TO BE CONSIDERED AS A POSSIBLE NATURAL 
FREQUENCY (RAD/SEC) 

THE WEIGHT OF MASS N (LBS) 


IMPLICIT REAL*8 ( A— H, Q-Z ) 
REAL*8 L 

CHARACTER*! KEY, YES 
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COMMON El (700) , EIO, FL (700) , L (700) , NPLACE, NRPM, NTOT, NW, RPM (10) 
1 , RRHUB, TITLE ( 15) , W ( 100) , WD, WF, WL, WT (700) , KEY 

DATA YES/’ Y’ / 

C 

C ASSIGN LOGICAL UNIT 1 TO ’VIBIN.DAT’ 

C 

0P£N(1, FILE=’ VIBIN. DAT’ , STATUS=’ OLD’ ) 

C 

C ASSIGN LOGICAL UNIT £ TO ’VIBOUT.DAT’ 

C 

OPEN (2, FILE=’ VIBOUT. DAT’ . STATUS=’ NEW’ ) 

C 

CJ READ DATA FROM VIBIN.DAT 

C 

READ (1, SO) TITLE 
90 FORMAT (6X, 15A4) 

READ (1,100) NTOT 
100 FORMAT (IX, 15) 

R£AD(1, 110) (L(N),£I (N> , WT (N) , N=1 , NTOT) 

1 10 FORMAT (IX, F10. 6, F20. 1, 7X, F15. 6) 

READ(1, 120) EIO 
120 FORMAT ( IX, F20. 1) 

READ (1,130) RRHUB 

130 FORMAT ( IX, F20. 3) 

C 

C READ ADDITIONAL DATA FROM TERMINAL KEYBOARD 

C 

131 WRIT£(*, 132) 

132 FORMAT (/, ’ TYPE 1 IF THE ANALYSIS IS FOR CHORDWISE’ , 

1 ’ VIBRATION OR’,/,’ TYPE £ IF THE ANALYSIS IS FOR’, 

2 ’ FOR FLAPWISE VIBRATION : ’ $) 

133 READ ( *, *>NA 

IF (NA. LT. 1. OR. NA. GT. 2) GO TO 131 
IF(NA. EQ. 1 ) GO TO 180 
135 WRITE (*, 140) 

140 FORMAT ( /, ’ TYPE IN THE TOTAL NUMBER OF ROTATIONAL ’, 

1 /, ’ SPEEDS TO BE CONSIDERED (MAX. 10) : ’ *> 

READ (*, * ) NRPM 

IF (NRPM. LT. 0. OR. NRPM. GT. 10) GO TO 135 
IF (NRPM. EQ. 0) GO TO 180 
C 

C ENTER ROTATION SPEEDS 

C 

DO 160 1 = 1, NRPM 
WRITE (*, 150) I 

150 FORMAT (/,’ TYPE IN ROTOR SPEED (REV/MIN) NO. ’,12,’ : ’ *) 

READ (*, *>RPM(I) 

160 CONTINUE 

GO TO 190 
C 

C NO ROTATIONS REQUESTED 

- C 

180 NRPM=1 

RPM ( 1 ) =0. 
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190 

200 

210 

220 

230 

240 

250 

C 

r*' 

u- 

c 

260 


WRITE (*, 200) 

FORMAT </,’ TYPE IN THE LOWER LIMIT OF THE FREQUENCY RANGE’,/, 

’ YOU WISH SEARCHED FOR NATURAL FREQUENCIES (RAD/SEC) : ’*) 
READ (*, *) WF 
WRITE <*, 210) 

FORMAT ( /, ’ TYPE IN THE UPPER LIMIT OF THE FREQUENCY RANGE’,/, 

’ YOU WISH SEARCHED FOR NATURAL FREQUENCIES (RAD/SEC) : ’*) 
READ <*, •*) WL 
WRITE (*, 220) 

FORMAT ( /, ’ TYPE IN THE FREQUENCY INCREMENT TO BE USED IN’,/, 

’ THE INITIAL SEARCH FOR NATURAL FREQUENCIES (RAD/SEC) : ’ *) 
READ (*, *) WD 
WRITE (*, 230) 

FORMAT < /, ’ TYPE IN THE NUMBER OF DIGITS OF ACCURACY TO THE’,/, 
’ RIGHT OF THE DECIMAL PGINT : ’ *) 

READ <*, * ) NPLACE 
WRITE <*, 240) 

FORMAT ( / /, ’ DO YOU WISH TO CALCULATE MODE SHAPES (Y/N) ? ’,*) 
READ (*, 250) KEY 
FORMAT ( A1 ) 

IF (KEY. NE. YES) GO TO 260 

ASSIGN LOGICAL UNIT 3 TO ’ MODES. DAT’ 

OPEN (3, FILE=’ MODES. DAT’ , STATUS=’ NEW’ ) 

CLOSE < 1 ) 

RETURN 

END 
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$nof loatcal Is A. 9 MODES - subroutine to calculate mode shapes 

^storage : £ 

SUBROUTINE MODES ( ISPEED) 


C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 


c: 

c 

c 

c 

c 

c 

c 

c 

c: 

c 

c 


c 

c 


c 

C 

c 

c 

c 

L 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


THIS SUBROUTINE CALCULATES THE MODE SHAPES If DESIRED. 


A ( I , J) THE MATRIX WHICH RELATES THE DEFLECTION, SLOPE, 

MOMENT AND SHEAR AT THE N-l MASS TO THOSE QUANTITIES 
AT MASS N 

ISPEED AN INDEX CORRESPONDING TO THE PARTICULAR ROTOR 
SPEED UNDER INVESTIGATION 

MO THE BENDING MOMENT AT THE FIXED END ASSUMING A UNIT 

DISPLACEMENT AT THE FREE END (INCH-LBS) 

NTOT THE TOTAL NUMBER OF MASSES 

NW THE NUMBER OF NATURAL FREQUENCIES IN THE RANGE FROM 

WF TO WL 

TITLE THE PROBLEM TITLE STATEMENT 

RPM(N) THE ROTOR ROTATIONAL SPEED (RPM) 

U(I,J) THE CONCATENATION OF CAD MATRICES RELATING DEFLECTION, 
SLOPE, MOMENT AND SHEAR AT THE FIXED END TO THOSE 
QUANTITIES AT MASS N 

UNEW ( I , J) SAME AS U(I,J> 

VO THE SHEAR FORCE AT THE FIXED END ASSUMING A UNIT 

DISPLACEMENT AT THE FREE END (LBS) 

WF THE FIRST VAuUE TO BE CONSIDERED AS A POSSIBLE NATURAL 

FREQUENCY (RAD/SEC) 

WL THE LAST VALUE TO BE CONSIDERED AS A POSSIBLE NATURAL 

FREQUENCY (RAD/SEC) 

W (NF) THE NATURAL FREQUENCIES OF VIBRATION (RAD/SEC) 

X (N) THE AXIAL LOCATION OF MASS N (INCHES) 

Y(N) THE NORMALIZED DISPLACEMENT AT MASS N. THAT IS, THE 
DISPLACEMENT AT MASS N ASSUMING A UNIT DISPLACEMENT 
AT THE FREE END OF THE BEAM 


4 
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c 

c 

c 

IMPLICIT REAL*8(A-H, 0-Z) 

REAL*8 L 

CHARACTER*! KEY, YES 

COMMON El (700) , EIO, FL (700) , L (700) , NPLACE. NRPM, NTOT, NW, RPM ( 10) 

1 , RRHUB, TITLE (15) , W(100) , WD, WF, WL, UT(700) , KEY 

REAL*8 A (4, 4) , U(4, 4) , UNEW(4, 4) , MO, Y (700) , X (700) , LTOT 
DATA YES/’Y’/ 

I F ( KEY. NE. YES ) RETURN 
IF ( ISPEED. NE. 1 ) GO TO 35 

10 FORMAT ( // * ** THE PROGRAM IS CALCULATING THE MODE SHAPES ** ’,//) 

C 

C WRITE THE TITLE, THE TOTAL NUMBER OF LUMPED MASSES 

C 1 AND THE AXIAL LOCATION OF EACH MASS. 

C 

WRITE (3, £0) 

£0 FORMAT(//, IX, ’ THE PROBLEM TITLE’ , /, IX, ’*****************’ ) 

WRITE <3, £5) TITLE 
£5 FORMAT ( /, 6X, 15A4) 

WRITE (3, £7) 

£7 FORMAT (//, IX, ’ THE TOTAL NUMBER OF MASSES’,/, 

1 IX, ’ **************************’ ) 

WRITE (3, 40) NTOT 
LTOT =0. 0 
DO 30 N=1 , NTOT 
LTOT=LTOT+L (N) 

30 X (N) =LTOT 

WRITE (3, 33) 

33 FORMAT (//, IX, ’ THE AXIAL LOCATION OF MASS N (INCHES)’,/, 

1 ’ ********************************#****’ , /) 

WRITE (3, 1E0) (X (N> , N=l, NTOT) 

C 

C WRITE THE TOTAL NUMBER OF MODES OF VIBRATION. 

C 

35 WRITE (3, 37) 

37 FORMAT (//, IX, ’ THE NUMBER OF NATURAL FREQUENCIES AT THE’ 

1 ,’ LISTED ROTOR SPEED (RPM)’,/, 

£ IX, ’ ***************************************************’ 

3 , ’ *************** ) 

WRITE (3, 40)NW, RPM(ISPEED) 

40 FORMAT (/, IX, 15, F10. 1 ) 

WRITE ( *, 10) 

C 

C CONSIDER ALL NATURAL FREQUENCIES 

C 

DO 130 NF=1,NW 
C 

C INITIALIZE CU3 (SET IT EQUAL TO THE IDENTITY MATRIX) 

C 

DO 50 1 = 1,4 
DO 50 J=l, 4 
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50 

C 

c 

c 

c 

c 

c 


60 
i 70 
C 
C 
C 

c 

c 


c 

c 

c 

c 

c 

c 


ao 

c 

c 

c 

c 

c 

c 

c 


so 

c 


U(I, J)=0. 

IF(I.EQ.J) U(I,J)=1. 

CONTINUE 

COMPUTE CUD = THE CONCATENATION OF THE CAD MATRICES FOR ALL MASSES 

DO 70 N=1 , NTOT 

FILL CAD FOR THE MASS N 

CALL FILL A (A, N, W (NF> ) 

CONCATENATE CAD WITH THE EXISTING CUD 
CALL MULT CA, 4, 4, 4, 4, U, 4, 4, 4, UNEW, 4, 4) 

SET CUD = CUNEWD 
DO 60 1 = 1,4 
DO 60 J=l, 4 
U(I, J)=UN£W(I, J) 

CONTINUE 

CALCULATE THE SHEAR AND MOMENT AT FIXED END ASSUMING A UNIT 
DISPLACEMENT AT THE FREE END (AUTOMATICALLY NORMALIZES THE MODE 
SHAPE WITH RESPECT TO THE DISPLACEMENT AT THE FREE END). 

M0=1. 0/ (U(l, 3) -U<4, 3>*U(1, 4) /U(4, 4) ) 

V0=— U<4, 3) *M0/U (4, 4) 

CALCULATE DISPLACEMENT AT MASS N IN TERMS OF THE MOMENT AND SHEAR AT 
THE FIXED END (MO AND VO RESPECTIVELY). 

REINITIALIZE CUD (SET IT EQUAL TO THE IDENTITY MATRIX) 

DO SO 1 = 1,4 
DO 80 J= 1,4 
U(I, J)=0. 

IF(I.EQ.J) U(I,J>=1. 

CONTINUE 

COMPUTE CUD = THE CONCATENATION OF THE CAD MATRICES FOR ALL 
MASSES UP TO AND INCLUDING MASS N 

DO 100 N= 1 , NTOT 

FILL CAD FOR THE MASS N 

CALL FILLA(A, N, W(NF) ) 

CONCATENATE CAD WITH THE EXISTING CUD 
CALL MULT (A, 4, 4, 4, 4, U, 4, 4, 4, UNEW, 4, 4) 

SET CUD = CUNEWD 
DO 90 1=1,4 
DO 90 J=1 , 4 
U ( I , J) =UNEW ( I, J ) 
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c 

c 

100 

105 


110 

115 


120 

130 


CALCULATE THE DISPLACEMENT OF MASS N 

V (N) =U (1,3) *M0+U (1,4) *V0 
WRITE (3, 105) 

FORMAT (//, IX, ’ THE NATURAL FREQUENCY OF VIBRATION (RAD/SEC)’,/, 

, IX, * ********************************************’ ) 

WRITE (3, 110) W <NF) 

FORMAT (/, 10X, F15. 5) 

WRITE (3, 115) 

FORMAT (//, IX, ’ THE NORMALIZED DISPLACEMENTS AT THE LUMPED MASSES’,/ 
IX, ’ *******#*#***#******#******************#*********’ , /) 

WRITE(3, 120) (Yd) , 1 = 1, NTOT) 

FORMAT (QF10. 5) 

CONTINUE 

RETURN , 

END 
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$nof loatcal Is A. 10 MULT -* subroutine to perform matrix multiplication 
$stc>rage:2 

SUBROUTINE MULT (A, ROWA, COLA, DROWA, DCOuA, B, COLB. DROWB, DCOLB, C, 

1 DROWC, DCOLC) 

C 

C THIS SUBROUTINE MULTIPLIES THE CAD MATRIX TIMES THE CBD 

C MATRIX AND STORES THE RESULTS IN THE CCD MATRIX. ALL 

C THREE MATRICES MUST BE UNIQUE <I.E. ONE CANNOT HAVE 


C 

p 

cad x 

CBD = 

CAl . 







L/ 

c 

p 

fid, j> 

THE 

FIRST MATRIX 





L/ 

c 

P 

B ( I , J ) 

THE 

SECOND 

MATRIX 





c 

r 

C ( I , J ) 

THE 

RESULT 

OF 

CAD X CBD = CCD 




c 

COLA 

THE 

NUMBER 

OF 

COLUMNS 

USED IN CAD 




c 


/ 

\ 








c 

p 

COLB 

THE 

NUMBER 

OF 

COLUMNS 

USED IN CBD 




w 

c 

DCOLA 

THE 

NUMBER 

OF 

COLUMNS 

DIMENSIONED FOR 

CAD 

IN THE 

CALLING 

c 

p 


PROGRAM 







c 

DCOLB 

THE 

NUMBER 

OF 

COLUMNS 

DIMENSIONED FOR 

CBD 

IN THE 

CALLING 

c 

p 


PROGRAM 







L 

c 

DCOLC 

THE 

NUMBER 

OF 

COLUMNS 

DIMENSIONED FOR 

CCD 

IN THE 

CALLING 

c 

p 


PROGRAM 







r*' 

L./ 

DROWA 

THE 

NUMBER 

OF 

ROWS DIMENSIONED FOR CAD 

IN 

THE CALLING 

c 

p 


PROGRAM 







L/ 

c 

DROWB 

THE 

NUMBER 

OF 

ROWS DIMENSIONED FOR CBD 

IN 

1 HE CALLING 

c 

p 


PROGRAM 







w 

c 

DROWC 

THE 

NUMBER 

OF 

ROWS DIMENSIONED FOR CCD 

IN 

THE CALLING 

c 

p 


PROGRAM 







L 

c 

ROWA 

THE 

NUMBER 

OF 

ROWS USED IN CAD 





C 

c 

REAL*Q A (DROWA, DCOLA) , B (DROWB, DCOLB) , C (DROWC, DCOLC) , VAL 
INTEGER COLA, COLB, DROWA, DCOLA, DROWB, DCOLB, DROWC, DCOLC, ROWA 
DO 10 1=1, ROWA 
DO 10 J= 1 , COLB 
VAL=0. 

DO Q K= 1 , COLA 

a VAL=VAL+A(I,K)*B(K, J> 

10 C ( I , J) =VAL 

RETURN 
END 
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$nof loatcal Is 


$storages£ 

SUBROUTINE NFREQ ( ISPEED) 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


THIS SUBROUTINE COMPUTES THE NATURAL FREQUENCIES IN THE RANGE 

FROM WF TO WL AND RETURNS THE RESULTS IN W(I> 

THE METHOD OF INTERVAL HALVING IS USED TO SEARCH FOR THE 

NATURAL FREQUENCIES. 

DET 1 THE DETERMINANT OF CUD AT THE STARTING POINT OF 
AN INTERVAL BOUNDING A NATURAL FREQUENCY 

DET£ THE DETERMINANT OF CUD AT THE ENDING POINT OF 
AN INTERVAL BOUNDING A NATURAL FREQUENCY 

DETM THE DETERMINANT OF CUD AT THE MIDPOINT OF 
THE INTERVAL BOUNDING A NATURAL FREQUENCY 

ISPEED AN INDEX CORRESPONDING TO THE PARTICULAR ROTOR 
SPEED UNDER INVESTIGATION 

NPLACE THE NUMBER OF DIGITS TO THE RIGHT OF THE DECIMAL POINT 

NW THE NUMBER OF NATURAL FREQUENCIES IN THE RANGE FROM 

WF TO WL 

W1 THE FREQUENCY AT THE BEGINNING OF AN INTERVAL OF LENGTH 

WD BOUNDING A NATURAL FREQUENCY (RAD/SEC) 

W1CUR THE FREQUENCY AT THE BEGINNING OF THE CURRENT INTERVAL 
OF LENGTH WDCUR BOUNDING A NATURAL FREQUENCY (RAD/SEC) 

W£CUR THE FREQUENCY AT THE END OF THE CURRENT INTERVAL OF 
LENGTH WDCUR BOUNDING A NATURAL FREQUENCY (RAD/SEC) 

W(I) THE NATURAL FREQUENCIES (RAD/SEC) 

WD THE INITIAL INTERVAL USED IN SEARCHING FOR NATURAL 

FREQUENCIES <RAD/SEC> 

WDCUR THE CURRENT FREQUENCY INTERVAL USED IN SEARCHING FOR 
A NATURAL FREQUENCY (RAD/SEC) 

WF THE FIRST VALUE TO BE CONSIDERED AS A POSSIBLE NATURAL 

FREQUENCY (RAD/SEC) 

WL THE LAST VALUE TO BE CONSIDERED AS A POSSIBLE NATURAL 

FREQUENCY < RAD/SEC) 

WMCUR THE FREQUENCY AT THE MIDPOINT OF THE INTERVAL 

BEING SEARCHED FOR A NATURAL FREQUENCY (RAD/SEC) 


IMPLICIT REAL*Q ( A— H, O— Z ) 
REALMS L 


126 



A. 11 NFREO - subroutine to determine natural frequencies 
CHARACTER*1 KEY 

COMMON El (700) , EIO, FL<700) , L<700) , NPLACE, NRPM, NTOT, NW, RPM(IO) 

1 , RRHUB, TITLE ( 15), W< 100), WD,WF,WL,WT (700), KEY 

WRITE (*, 3) 

3 FORMAT ( / //, £X, * ** THE PROGRAM IS NOW DETERMINING THE NATURAL 

1 /, 5X, ’ FREQUENCIES OF VIBRATION ** ’ ) 

WRITE (*, 4) RPM ( ISPEED) 

4 FORMAT ( / // ’ NATURAL FREQUENCIES AT A ROTOR SPEED OF ’ , F6. 1,, 

1 ’ RPM » ) 

WRITE (*, 5) 

5 FORMAT ( ’ ************************************************** » ) 
C 

C CALCULATE FORCES DUE TO ROTATION 

C 

CALL FORCE (ISPEED) 

C 

C INITIALIZE VARIABLES 

C 

NW=0 

W1=WF 

6 W1CUR=W1 
DET1=DET (W1 ) 

10 W2CUR=W1CUR+WD 

C 

C CHECK IF DONE 

C 

IF (W1CUR. GE. WL) RETURN 
IF (W2CUR. GT. WL) W2CUR=WL 
C 

C NOT DONE, CONTINUE SEARCHING 

C 

DET2=DET (W2CUR) 

C 

C CHECK IF A NATURAL FREQUENCY OCCURS BETWEEN W1CUR AND W2CUR, 

C IS AT W1CUR OR W2CUR, OR IS NOT IN THE INTERVAL W1CUR TO W2CUR 

C 

15 IF (DET1*DET2) 21,30,40 

C 

C A NATURAL FREQUENCY OCCURS BETWEEN W1CUR AND W2CUR 

C 

C COMPUTE THE SIZE OF INTERVAL IN ORDER 

C TO CHECK IF DESIRED ACCURACY HAS BEEN REACHED 

C 

21 WDCUR=W2CUR— W 1 CUR 

IF ( WDCUR. LT. 10. ** (-NPlACE-1 ) ) GO TO 25 
C 

C WDCUR TOO LARGE, CONSIDER FREQUENCY AT INTERVAL MIDPOINT 

C 

WMCUR= (W1CUR+W2CUR) /2. 

DETM=DET (WMCUR) 

C 

C DECIDE WHICH ENDPOINT, W1CUR OR W2CUR, 

C TO REDEFINE 
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c 

c 

c 


c 


es 


c 

c 

c 

c 

£3 


C 

C 

C 


1000 


c 

c 

c 

30 


C 

C 

C 


c 

c 

c 


c 

c 

c 

35 


c 

c 


c 


c 

c 

c 

40 


IF(DETM*D£T1) 2£, £3, £3 

DETM IS ON OPPOSITE SIDE OF AXIS FROM DET 1 SO REDEFINE W£CUR 

DET£=DETM 
W£CUR=WMCUR 
GO TO 15 

DETM IS ON THE SAME SIDE OF AXIS AS DET1 OR IS ZERO 
SO REDFINE W1CUR 

DET1=DETM 
W1CUR=WMCUR 
GO TO 15 

FOUND ROOT 

NW=NW+1 

W(NW)=(W1CUR+W£CUR) /£. 

WRITE (*, 1000) W <NW) 

FORMAT ( / ’ W = * , F 1 5. 3, ’ RAD/SEC ’) 

Wl=Wl+WD 
GO TO & 

ONE OF THE ENDPOINTS WAS A NATURAL FREQUENCY 

IF <DET£. EQ. O. > GO TO 35 

W1CUR IS A NATURAL FREQUENCY, STORE IT 

NW=NW+1 

W(NW)=W1CUR 

WRITE <*, 1000) W(NW) 

MOVE TO NEXT INTERVAL 

W1=W1+WD 
GO TO 6 

W£C'UR MUST BE A NATURAL FREQUENCY, STORE IT 

NW=NW+1 

W(NW)=W£CUR 

WRITE (*, 1000) W(NW) 

MOVE TO NEXT INTERVAL 

Wl=Wl+£, *WD 
GO TO & 

NO NATURAL FREQUENCY FOUND, MOVE TO NEXT RANGE 

W1CUR=W£CUR 
W1=W1CUR 
DET1=DET2 
GO TO 10 
END 
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A. 12 OUTPUT - subroutine to perform all output functions 

$nof loatcal Is 
$stor'ages£ 

SUBROUTINE OUTPUT ( ISPEED) 

C 

C THIS SUBROUTINE OUTPUTS THE INPUT DATA AND NATURAL FREQUENCIES 

C TO FILE ’VIBOUT.DAT’ 

C 

C 

C El <N> THE PRODUCT OF MODULUS OF ELASTICITY AND AREA MOMENT 

C OF INERTIA AT MASS N (LBS* INCH* INCH) 

C 

C EIO THE PRODUCT OF MODULUS OF ELASTICITY AND AREA MOMENT 

C OF INERTIA AT ZERO OR THE FIXED END (LBS*INCH*INCH) 

C 

C ISPEED AN INDEX CORRESPONDING TO THE PARTICULAR ROTOR 

C SPEED UNDER INVESTIGATION 

C 

C L <N) THE LENGTH OF SECTION TO LEFT OF MASS N (INCHES) 

C 

C NPLACE THE NUMBER OF DIGITS REQUIRED TO THE RIGHT OF THE 

C DECIMAL POINT 

C 

C NTOT THE TOTAL NUMBER MASSES 

C 

C NW THE NUMBER OF NATURAL FREQUENCIES IN THE RANGE FROM 

C WD TO WL 

C 

C RPM(N) THE ROTOR ROTATIONAL SPEED ( REV/MIN) 

C 

C RRHUB THE RADIUS OF THE RIGID HUB (INCHES) 

C 

C TITLE THE PROBLEM TITLE STATEMENT 

C 

C W(I> THE NATURAL FREQUENCIES OF VIBRATION (RAD/SEC) 

C 

C WD THE INITIAL DELTA USED IN SEARCHING FOR NATURAL 

C FREQUENCIES (RAD/SEC) 

C 

C WF THE FIRST VALUE TO BE CONSIDERED AS A POSSIBLE NATURAL 

C FREQUENCY (RAD/SEC) 

C 

C WL THE LAST VALUE TO BE CONSIDERED AS A POSSIBLE NATURAL 

C FREQUENCY (RAD/SEC) 

C 

C WT (N) THE WEIGHT OF MASS N (LBS) 

C 

C 

C 

C 

IMPLICIT REAL*8 ( A— H, O— Z ) 

REAL*8 L 
CHARACTER*! KEY 
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c 

c 

c 

99 

100 
110 


120 


140 

150 

160 


165 

170 


175 

180 


300 

301 

310 

31 1 

320 

321 

330 

331 

340 

341 
190 


COMMON El (700) , EIO, FL(700) , L(700) , NPLACE, NRPM, NTDT, NW, RPM(IO) 

, RRHUB, TITLE ( 15) , W (100) , WD, WF, WL, WT (700) , KEY 
IF ( ISPEED. NE. 1 ) GO TO 175 

WRITE OUT DATA READ FROM FILE ’VIBIN.DAT’ 

WRITE <2, 99) TITLE 
FORMAT <6X, 15A4) 

WRITE (2, 100) NTOT 

FORMAT (//’ THE NUMBER OF LUMPED MASSES = ’,15) 

WRITE (2, 110) 

FORMAT ( //’ SECTION LENGTH SECTION RIGIDITY’, 

’ LUMPED WEIGHT ’ ) 

WRITE (2, 120) 

FORMAT </’ (INCH) (LBS-IN*IN) ’ , 

’ (LBS) ’ ) 

DO 140 N=1 , NTOT 
WRITE (2, 150) L(N),EI(N),WT(N> 

FORMAT (/, IX, F10. 6, 3X, F20. 1,3X,F15. 6) 

WRITE (2, 160) EIO 

FORMAT (//,’ THE SECTION RIGIDITY (AT X=0. O) = ’ , F20. 1, 

’ (LBS— IN*IN) ’ ) 

WRITE (2, 165) RRHUB 

FORMAT <//,'* THE RADIUS OF THE RIGID HUB = ’,F10. 3, ’ (IN) ’) 
WRITE(2, 170) WF, WL, WD 
FORMAT (//,’ THE FREQUENCY RANGE SEARCHED FOR NATURAL FREQUENCIES ’,/ 
’ STARTING AT:’, FIS. 5,’ AND ENDING AT: ’ , F15. 5, ’ (RAD/SEC)’, 

//, ’ THE INITIAL FREQUENCY INCREMENT USED IN THE SEARCH = »,F15. 5 
,’ (RAD/SEC)’) 

FACTOR* 1 . / (2. *3. 14159265) 

WRITE (2, 180) RPM( ISPEED) 

FORMAT (////, 8X, ’ NATURAL FREQUENCIES AT A ROTOR SPEED OF ’ , F6. 1, 

’ RPM’ 

, /, 8X, ’ *************"*************************************’ , //, 

6X, ’ (RAD/SEC) ’ , 16X, ’ (HZ) ’ , //) 

II=NPLACE+1 
IF (NPLACE. LT. O) 11=1 
IF (NPuACE. GT. 4) 1 1=5 
DO 190 1=1, NW 

GO TO <300, 310, 320, 330, 340) II 
WRITE (2, 301) I, W(I) , W(I)*FACTOR 
FORMAT (/, IX, 14, F10. O, 10X, F10. O) 

GO TO 190 

WRITE (2, 311) I , W < I ) , W < I ) *FACTOR 
FORMAT (/, IX, 14, F10. 1, 10X, F10. 1 > 

GO TO 190 

WRITE (2, 321) I,W(I),W(I) *F ACTOR 
FORMAT (/, IX, 14, F10. 2, 10X, F10. 2) 

GO TO 190 

WRITE (2, 331) I, W(I) , W(I)*FACTOR 
FORMAT (/, IX, 14, F10. 3, 10X, F10. 3) 

GO TO 190 

WRITE (2, 341) I , W ( I ) , W < I ) *FACTOR 
FORMAT ( /, IX, 14, F10. 4, 10X, F10. 4) 

CONTINUE 

RETURN 

END 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


50 

c 

c 

c 

c 

c 

c 

c 

c 

c 

100 

c 

c 

c 

110 

120 


c 

c 

c 

130 

c 

c 

c 

140 

145 


A. 13 MODPLT - mode shape plotting code 

THIS PROGRAM Pi_OTS THE MODE SHAPES FOR THE LUMPED MASS MODEL 
OF THE BEAM VIBRATION 


AMSF 

ITERM 

NC 

NTOT 

NW 

TITLE 

W 

X (NTOT) 
Y (NTOT) 


THE AMPLITUDE SCALE FACTOR 

THE TERMINAL NUMBER (4012 OR 4014) FOR GRAPHICS 

THE NUMBER OF CHARACTERS IN THE TITLE 

THE TOTAL NUMBER OF LUMPED MASSES 

THE TOTAL NUMBER OF NATURAL FREQUENCIES 

THE TITLE TO BE PRINTED AT THE TOP OF THE PLOT 

THE CURRENT NATURAL FREQUENCY 

THE X COORDINATES OF THE LUMPED MASSES 

THE DEFLECTIONS OF ThE LUMPED MASSES 


BYTE TITLE 

COMMON AMSF, NC, NTOT, NW, TITLE (60) , W, X (900) , Y(900) 

DATA AMSF/ 10. / 

OPEN (FILE-* MODES. DAT’, UNIT=1 , TYP£=’ NEW’ ) 

CALL START 
WRITE (*, 50) 

FORMAT ( ’ ENTER TYPE OF DISPLAY TERMINAL (4012 OR 4014) ’) 
READ (*, *) ITERM 


READ IN NUMBER OF MASSES, NTOT, AND THE NUMBER OF MODES OF 
VIBRATION 


CALL MINPUT ( 1 ) 

READ NATURAL FREQUENCY AND INITIAL MODE SHAPE (IT MAY NOT BE 
FIRST MODE IF A NON-ZERO INITIAL FREQUENCY WAS GIVEN TO THE 
BEAM PROGRAM). 

CALL MINPUT (2) 

SCALE DRAWING 

WRITE(*, 120) AMSF 
FORMAT (’ 1 = DRAW SKETCH’ 

,/, ’ 2 = CHANGE AMPLITUDE SCALE FACTOR, CURRENTLY = ’ , F6. 2 

,/, ’ 3 = DRAW FINISHED DRAWING (S) ’> 

READ (*, *> I 

IF ( I. LT. 1. OR. I. GT. 3) GO TO 110 
GO TO (130, 140, 150) , I 

DRAW SKETCH 


CALL DRAW (ITERM, 1) 

GO TO HO 

CHANGE SCALE FACTOR 
WRITE (*, 145) 

FORMAT ( ’ ENTER NEW SCALE FACTOR TO BE APPLIED TO AMPLITUDE ’) 
READ (*, *) AMSF 


131 



go to no 
c 

C DRAW FINISHED DRAWING (S) FOR ALL FREQUENCIES USING SAME SCALE 

C FACTOR 

C 

ISO IF (NW. LE. 0) STOP 

DO 160 1=1, NW 
CALL DRAW ( ITERM, £> 

C READ IN NEXT MODE SHAPE IF THERE IS ANOTHER MODE SHAPE IN DATA 

C FILE 

I F < 1 - LT . NW) CALL MINPUT<£> 

160 CONTINUE 

C 

C STOP 

C 

STOP 

END 
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A. 14 DRAW - subroutine to draw mode shapes 


SUBROUTINE DRAW ( ITERM, ITYPE) 

C 

C THIS SUBROUTINE DRAWS ft SIMPLE SKETCH OF THE MODE SHAPE 

C (ITYPE = 1) FULL DRftWING COMPLETE WITH TITLE BLOCK (ITYPE = £) . 

C flMSF THE AMPLITUDE SCftLE FACTOR 

C CHARX THE X SIZE OF THE INDIVIDUAL CHARACTERS TO BE DISPLAYED 

C DAT THE CHARACTERS FOR THE CURRENT DATE 

C LFTMAR THE SIZE OF THE LEFT MARGIN ON THE SCREEN 

C MODES THE CURRENT MODE NUMBER 

C NC THE NUMBER OF CHARACTERS IN THE TITLE 

C NTOT THE TOTAL NUMBER OF LUMPED MASSES 

C NW THE TOTAL NUMBER OF NATURAL FREQUENCIES 

C RGTMAR THE SIZE OF THE RIGHT MARGIN ON THE SCREEN 

C SIGN THE SIGN OF THE AMPLITUDE AT THE PREVIOUS LUMPED MASS 

C TITLE THE TITLE TO BE PRINTED 

C W THE CURRENT NATURAL FREQUENCY 

C Y (NTOT) THE X COORDINATES OF THE LUMPED MASSES 

C XNODES THE X COORDINATES OF THE NODES 

C XORG THE X SCREEN COORDINATE FOR THE AXES ORGIN 

C XSCALE THE SCALE FACTOR TO BE APPLIED TO ALL X COORDINATES 

C XSCRIM THE WIDTH OF THE SCREEN IN INCHES 

C XSTART THE X SCREEN COORDINATE WHERE TITLE BLOCK BEGINS 

C Y (NTOT) THE Y DEFELECTIONS OF THE LUMPED MASSES 

C YMAX THE MAXIMUM DEFLECTION OF THE LUMPED MASSES 

C YMIN THE MINIMUM DEFLECTION OF THE LUMPED MASSES 

C YORG THE Y SCREEN COORDINATE FOR THE AXES ORIGIN 

C YP THE SCALED Y DEFLECTION 

C YPOLD THE PREVIOUS SCALED Y DEFLECTION 

C YSCRN THE HEIGHT OT THE SCREEN IN INCHES 

C 

BYTE TITLE, DAT (S) 

COMMON AMSF, NC, NTOT, NW, TITLE (SO) , W, X (900) , Y(900) 

REAL LFTMAR, XNODES (100) 

DATA LFTMAR, RGTMAR/. 5, . 5/ 

WRITE (*,100) 

100 FORMAT ( ’ PUSH THE RETURN KEY TO DISPLAY NEXT PlGT’ , /, 

1 » THEN AFTER PLOT IS DISPLAYED PUSH RETURN KEY FOR NEXT PLOT’ ) 

READ (*, *) I 
CALL PLOTS (I TERM) 

CALL ERASE 

SET TERMINAL VALUES 

IF ( I TERM. EQ. 4014) GO TO 120 
C 4012 

XSCRN=7. 307 
YSCRN=5. 564 
CHARX=. 1085 
GO TO 130 
C 4014 

ISO XSCRN=14. 62 

YSCRN=1 1 . 14 
CHARX =0. 196 
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c 

C SCALE X 

C 

130 XSCALE= (XSCRN— LFTMAR— RGTMAR) /X (MTOT) 

C 

C LOCATE AXIS 

C 

YORG= (YSCRN-. 5-1. O) /£. +1. 

XORG=LFTMAR 

C 

C PLOT SOLID LINE MODE SHAPE, LOOK FOR AXIS CROSSINGS TO DETERMINE 

C MODE NUMBER AND LOOK FOR MAX AND MiN AMPLITUDE VALUES 

C 

M0DES=1 

SIGN=0 

CALL PLOT (XORG, YORG, -3) 

YMAX=0 

YMIN=0 

DO £00 N=1 , NTOT 
YP=Y (N) *AMSF 
YMAX=AMAX1 (YP, YMAX) 

YMIN=AMI!M1 (YP, YMIN) 

IF (SIGN. NE. O) GO TO 180 

SIGN=YP 

GO TO 190 

C CHECK FOR CROSSING 

180 IF (YP*SIGN) 185,185,190 

C CROSSING FOUND 

185 M0DES=M0DES+1 

C COMPUTE NODE X COORDINATE LOCATION 

XNODES (MODES) =X (N-l > + (YPOLD-YP) / (YPGLD*(X (N) —X (N— 1 ) ) ) 

SIGN=YP 

C NO CROSSING FOUND, PLOT POINT 

190 CALL PLOT ( X (N) *X SCALE, YP, £) 

YPOLD=YP 
£00 CONTINUE 

C 

C PLOT AXIS 

C 

CALL PLOT (X (NTOT) *XSCALE,0. ,3) 

CALL PLOT (0. , 0. , £) 

CALL PLOT (O. , YMAX, £> 

CALL PLOT (O. , YMIN, £) 

C CHECK IF FULL DRAWING IS TO BE DONE 

IF ( ITYPE. EQ. £) GO TO 300 
C NO, ONLY SKETCH 

CALL PLOT (XSCRN— XORG, YSCRN— YORG, 999) 

READ (*, *) I 
CALL ERASE 
RETURN 
C 

C DRAW FULL DRAWING 

C 

C BEGIN BY DRAWING DASHED ENVELOPE (SOLID LINE DRAWN ABOVE) 

300 IPEN=£ 
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CALL PLOT <0. , 0. , 3) 

DO 320 N=l, NTOT 

CALL Pl.OT ( X (N) *X SCALE, -Y(N)*AMSF, IPEN) 

C SWITCH FROM SOLID TO DRAWN LINE FOR DASHED EFFECT 

IF < IPEN. EQ. 2) 60 TO 310 
I PEN=2 
GO TO 320 
310 I PEN=3 

320 CONTINUE 

C 

C PLOT X/L LOCATIONS OF NODES 

C 

C PLOT HEADER 

CALL PLOT (-4. *CHARX, YMIN-. 2, 3) 

CALL STRINGS X/L:’ , 4, 4) 

DO 330 1=1, MODES 
XX=XNODES ( I ) *XSCALE 
CALL PLOT (XX, —0. 05, 3> 

CALL PLOT (XX, YMIN, 2) 

CALL NUMBER (XX-2. 5*CHARX, YMIN-. 2, O. 11, XNODES < I ) / X < NTOT ) , O. , 
1 ’ (FE. 3) ’ ) 

330 CONTINUE 

PLOT TITLE AND TITLE BLOCK 

RESET ORIGIN 

CALL PLOT (-XORG, -YORG, —3) 

C PLOT TITLE 

CALL PLOT (XSCRN-FLOAT (NC+l ) /2. *CHARX, YSCRN-. 5, 3) 

CALL STRING (TITLE, NC, 4) 

C PLOT TITLE BLOCK 

XSTART= (XSCRN-4. 5) /2. 

CALL PLOT(XSTART, 0. , 3) 

CALL PLOT (XSTART+5. , O. , £> 

CALL PLOT (XSTART+5. , 1. , 2) 

CALL PLOT (XSTART, 1. , 2) 

CALL PLOT ( XSTART, O. ,2) 

CALL PLOT (XSTART, . 3, 3) 

CALL PLOT ( XSTART+5. ,. 3, 2) 

CALL PLOT (XSTART+5. , 1. , 3) 

CALL PLOT (XSTART+5. , 1. , 3) 

CALL PLOT (XSTART+5. , 0. , 2) 

CALL PLOT < XSTART+1 . 75, 1 . , 3) 

CALL PLOT (XSTART+1. 75, 0. , £> 

CALL PlOT (XSTART+3. 25, 1. , 3) 

CALL PLOT (XSTART+3. 25, 0. , 2) 

CALL PLOT (XSTART+. 1, .6,3) 

CALL STRING (’MODE’ ,4, 1) 

CALL PLOT ( XSTART+. 6, . 6, 3) 

CALL STRING (’NAT. FREQ.’, 10,1) 

CALL PLOT (XSTART+. 6, . 3, 3) 

CALL STRING (’ (RAD/SEC) 9, 1 ) 

CALL PlOT(XSTART+1. 85, . 6, 3) 

CALL STIRNG (’ NO. OF MASSES’ , 13, 5 ) 

CALL Pt_OT (XSTART+3. 35, . 6, 3) 
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CALL STRING (’ DATE’ ,7,1) 

CALL NUMBER ( XSTART+. 1, . 1,0. 1 1 , MODES, 0. , ’ (13) » ) 
CALL NUMBER (XSTART+. 6, . 1, 0. 11, W, 0. , ’ (F10. 3) ’ ) 
COLL NUMBER (XSTART+l. B5, .1,0. 11, NTOT, O. , » (17) ’ ) 
CALL DATE (DAT) 

CALL PLOT (XSTART+3. 35, . 1, 3) 

CALL STRING (DAT, 9, 1 ) 

CALL PLOT (XSCRN, YSCRN, 999) 

READ ( *, * ) I 
CALL ERASE 
RETURN 
END 
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c 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 
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A. 15 MINPUT «* subroutine to perform all input functions 
SUBROUTINE MINPUT ( ITYPE) 

THIS SUBROUTINE READS DATA FROM DATA FILE ’ MODES. DAT’ . 


NC 

NTOT 

NW 

TITLE 

W 

X (NTOT) 
Y (NTOT) 


THE NUMBER OF CHARACTERS IN TITLE 
THE TOTAL NUMBER OF MASSES 

THE TOTAL NUMBE OF NATURAL FREQUENCIES 

THE TITLE TO BE PRINTED 

THE CURRENT NATURAL FREQUENCY BEING CONSIDERED 

THE X COORDINATES OF THE LUMPED MASSES 
THE Y DEFLECTIONS AT THE LUMPED MASSES 


IF ITYPE = 1 THE TITLE, NUMBER QF LUMPED MASSES, THEIR 
X COORDINATES AND THE NUMBER OF NATURAL FREQUENCIES COMPUTED 
ARE READ. 

IF ITYPE = £, THE NATURAL FREQUENCY AND MODE SHAPE AMPLITUDES 
AT THE LUMPED MASSES ARE READ. 

BYTE TITLE 

COMMON AMSF, NC, NTOT, NW, TITLE (60) , W, X (300) , Y (900) 

BRANCH ON TYPE OF READ 


GO TO (100, £00) , ITYPE 
C 

C READ TITLE, NUMBER OF LUMPED MASSES, THEIR X LOCATIONS AND THE 

C NUMBER OF MODES ARE READ. 

C 

100 READ ( 1 , 110) NC, (TITLE(I) , 1=1, NC) 

110 FORMAT (6X, Q., 60A1) 

READ ( 1, 1£0) NTOT 
ISO FORMAT (IX, 15) 

READ (1,130) (X (N> , N=l, NTOT) 

130 FORMAT (8F 10. 5) 

READ ( 1 , 1£0) NW 
RETURN 
C 

C READ NATURAL FREQUENCY AND MODE SHAPE 

C 

£00 READ ( 1 , £10) W 

£10 FORMAT ( 10X, F15. 5) 

READ (1, 130) (Y(N) , N=l, NTOT) 

RETURN 

END 


«U.S. GOVERNMENT PRINTING OFFICE: 1986-646-114-20432 
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